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This  report  considers  techniques  used  to  solve  the  Ravier  field  equations  of 
viscoelasticity  in  the  Kelvin-Voigt  or  Maxvell  models,  for  cylindrical  or  spher 
ical  geometries.  Introducing  scalar  and  vector  potentials  into  the  viscoelas- 
ticity equations  formulation,  ultimately  yields  **telegraph-type’wl3partial 
differential  equations  governing  those  potentials.  For  harmonic  time-dependencel, 
these  reduce  to  scalar  and  vector  Helmholtz 1 s equations  with  complex  propagation 
constants.  These  constants  are  shovn  tc  be  related  to  the  viscoelastic  material 
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constants  in  a more  or  less  complicated  fashion  depending  on  the  viscoelastic 
model  used.  The  stresses,  strains  and  displacements  are  then  found  from  thesi 
potentials  for  a dozen  cases  of  interest  in  those  two  coordinate  systems.  Th< 
formulation  resembles  that  of  electrodynamics  in  a Coulomb  gauge 

' Jhe  above  information  is  vital  to  set-up  and  solve  various  kinds  of 
boundary-value-problems  of  dynamic  viscoelasticity  which  appear  when  studying 
cases  of  acoustic  scattering  from  sound-absorbing  structures ^problems  we  are 
now  addressing.  The  analysis  is  summarized  in  two  large  TablgP£t-§£t-up  in  a 
conveniently  accessible  form.  Remarks  on  "complex-moduli"  af’d  examined  in  a 
final  section  under  the  light  of  the  viscoelasticity  "Correspondence  Theorem" 
and  a list  of  recommendations  and  conclusions  is  given  at  the  aid. 
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METHODS  FOR  SOLVING  THE  VISCOELASTICITY  EQUATIONS  FOR  CYLINDER 
AND  SPHERE  PROBLEMS 


This  report  describes  techniques  used  to  solve  the  field  equations  of 
dynamic  viscoelasticity  in  tvo  geometries,  tvo  viscoelastic  models,  and  a dozen 
different  situations  of  interest.  This  information  is  needed  to  study  acoustic 
scattering  from  sound-absorbing  structures  of  cylindrical  and  spherical  shapes, 
which  are  cases  presently  under  study  by  the  author.  This  report  sets  up  most 
of  the  theoretical  viscoelasticity  foundations  needed  to  deal  with  the  other 
sound  scattering  problems  under  study. 

This  work  is  continuing  and  it  was  done  as  part  of  an  NSWC  project  entitled 
''Acoustical  Properties  of  Ordnance  Materials",  Task  No.  MAT-Q3L-000/ZROO-001-Q10 
Problem  127,  which  deals  with  acoustic  scattering  from  objects  covered  with 
viscoelastic  materials.  This  is  a progress  report  describing  work  done  during 
FY  76. 
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I.  INTRODUCTION  The  basic  viscoelastic  models. 

The  deformation  of  a viscoelastic  solid  under  any  kind  of  external  loads  is 
usually  studied  by  means  of  viscoelastic  models.  Two  basic  models  commonly  used 
are  associated  with  the  names  of  Kelvin1 -Voigt 2 and  Maxwell3.  In  the  Kelvin- 
Voigt  model,  the  elastic  and  viscous  properties  of  each  material  point  (or  par- 
ticle) of  the  body,  which  can  be  respectively  represented  by  a spring  and  a dash- 
pot,  are  assumed  connected  in  parallel.  In  the  Maxwell  model  they  are  assumed 
connected  in  series.  We  rejeat  that  this  description  applies  at  each  point  in 
the  viscoelastic  solid  and  it  is  as  if  the  body  contained  a continuous  distri- 
bution of  damped  oscillators,  (viz,  elementary  mass-spring-dashpot  systems)  one 
at  each  one  of  its  material  points.  As  if  it  were  not  clear  enough  already,  the 
entire  "Kelvin-Voigt  solid’  can  not  be  replaced  by  ONE  spring  connected  in  paral- 
lel with  ONE  dashpot.  This  simplistic  view  of  a deformable  solid  as  a single 
particle,  may  be  useful  in  some  other  elementary  context  such  as  that  used  when 
one  treats  a body  as  a particle,  but  it  is  of  no  use  in  viscoelasticity.  Other- 
wise viscoelasticity  could  not  be  viewed  as  a field-theory  capable  of  describing 
stress  and  displacement  fields  at  each  point  in  a body,  since  by  that  oversimpli- 
fication, the  body  has  been  reduced  to  a particle.  We  emphasize  this  rather 
trivial,  but  quickly  forgotten  point,  because  it  is  common  to  find  authors  who 
try  to  use,  say,  the  "Kelvin  model"  as  a single  spring  in  parallel  with  a single 
dashpot,  only  to  find  that  the  model  is  "no-good"  and  that  they  must  go  to  "more- 
degrees-of-freedon  systems",  such  as  three  "Kelvin-model 5"  in  series,  or  other 
similar  configurations,  to  obtain  meaningful  results.  It  is  clear  that  this  ap- 
proach does  not  give  the  Kelvin  model,  as  it  truly  is,  even  a char.ce  to  H/ork". 
These  authors  have  replaced  the  continuous  field-equations  of  viscoelasticity **, 
by  a set  of  three  ordinary  second-order  differential  equations  of  the  simple 


1 William  Thomson  (Lord  Kelvin) , b.  at  Belfa^-,  .bu  ; d.  near  Glasgow,  1907. 
Professor  at  Glascow  Univ,  16U6— 1809 • Buried  at  Westminster  Abbey,  London 
(near  Newton's  tomb). 

1 Voldemar  Voigt . b.  at  Leipzig  1850;  d.  at  Gottingen  1919*  Professor  of 
Mechanics  at  Gottingen  Univ,  Germany. 

1 James  C.  Maxwell,  b.  Edinburgh,  1@31;  d.  Cambridge,  l8?9*  Professor  of 
Physics  at  King's  College  and  at  Cambridge  Univ,  1860-1879.  Founder  of  the 
Cavendish  Lab  at  Cambridge. 

**  ie,  a set  of  three  scalar  partial  differential  equations  hopelessly  coupled 
and  non-linear,  governing  the  displacement-field  in  the  body,  which  can  be 
linearised  for  small  deformations  and  small  deformation-gradients,  and  the 
"linear"  theory  of  viscoelasticity  then  results.  These  equations  are  called 
the  Navier  equations  of  viscoelasti 'ity. 
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damped-oscillator  type.  This  replacement  is,  in  no  way,  equivalent  to  solving  the 
field-equations  of  Navier. 


The  basic  point  of  the  continuum  field-theory  approach  to  viscoelasticity  is 
that  the  model  assumptions,  that  the  spring  and  dashpot  are  connected  in  series 
or  in  parallel  at  each  point  in  the  body,  are  immediately  reflected  in  the  fact 
that  the  field-equations  resulting  from  either  one  of  those  models  turn  out  to  be 
different.  To  fix  the  ideas  we  now  give  the  linearized  form  of  the  field-equations 
for  both  these  models. 


1)  Kelvln-Voigt  Model  (parallel): 


i] 


• z>  - 4 


»fs 

at2 


il)  Maxwell  Model  (series): 


2+ 


+ U(f  .5)  i i_ a 


at* 


3a 

26 

y. 


72u  + 


3*u  u.  3a  + 26  3u  .„2  „ 

7?  ~7~~  3t  8 p 


^ Here  0 is  the  material  density,  X , u are  the  elastic  Lame  constants, 
u is  the  vector-displacement  field,  and  eX  ,ey  (or  a,  6 in  the.Maxvell  model) 
are  the  viscosity  coefficients.  As  given  a&jve'f  in  differential  operator  form, 
these  equations  hold  in  any  coordinate  system.  These  are  the  Havier  equations  one 
must  solve  in  the  viscoelastic  body.  It  is  possible  to  derive  eqs.  (i)  and  (ii) 
starting  from  the  basic  idea  that  the  spring  and  dashpot  at  each  material  point 
are  connected  either  in  parallel  or  in  series  respectively.  Once  the  displacement 
field  components  are  found  by  solving  (i)  or  (ii)  with  suitable  boundary  conditions, 
one  can  then  find  the  stresses  from  them. 


In  the  absence  of  viscosity  (ie,  Xy  «*  0,  u„  ® 0 for  the  Kelvin  model  or 
a * 0,  8 «*  0 for  the  Maxwell  model)  both  field  equations  (i)  and  (ii)  given  above 
reduce  to, 


iii)  The  field  equations  of  linear  dynamic  elasticity: 


» 
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since  now  there  are  only  springs  and  no  dashpots  at  each  material  point.  The 
scalar  components  of  this  vector  equation  are  the  Navier  equations  of  elasticity 
in  the  absence  of  body- forces  as  found  by  Navier5  for  the  elastic  body.  Equations 
(i)  and  (ii)  are  also  called  the  Navier  equations  (of  viscoelasticity)  by  exten- 
sion, since  Navier  never  worked  with  viscoelastic  solids. 

In  most  materials  of  interest,  the  "constants"  X^,  and  the  "coefficients" 
X , y are  really  not  constants,  but  frequency-dependent  parameters.  Thus,  few 
materials  are  completely  describable  by  the  Kelvin  or  the  Maxwell  models  in  the 
continuum  sense  of  eqs.  (i)  and  (ii).  The  material  behavior  of  viscoelastic  sub- 
stances tends  to  follow  one  or  the  other  model  in  different  regions  of  the  para- 
meters involved,  say,  frequency  among  others.  This  means  that  in  general,  these 
models  of  field-equations  (i)  and  (ii)  are  quite  "good".  Rubbers  at  low  frequen- 
cies are  known  to  be  well  described  by  the  Kelvin  model  (i).  Pulse-tube  measure- 
ments exploit  this  factual  observation.  Since  these  models  are  not  perfect,  re- 
searchers in  this  field  have  proposed  more  complicated  models. 


It  is  not  hard  to  see  that  various  Maxwell  elements  in  series  at  each  point 
in  the  body  have  the  properties  of  a single  Maxwell  element  with  equivalent  spring 


and  dashpot  constants  given  by 


and  1/n 

eq 


respectively.  Various  Kelvin  elements  in  parallel  at  each  point  in  the  body  have 
the  properties  of  a single  Kelvin  element  vich  ^ aQl*  ''eq  " "i* 


On  the  other  hand,  Kelvin  elements  in  series,  or  Maxwell  elements  in  parallel, 
have  more  complicated  properties.  In  their  desire  to  generalise  the  basic  models 
(i)  and  (ii),  researchers  have  invented  the  so  called  "standard  viscoelastic 
model".  It  consists  of  a Maxwell  element  in  parallel  with  a Kelvin  element  at 
each  point  in  the  body.  The  linearised  field  equations  which  result  in  this  model 
when  two  springs  and  two  dashpots  are  connected  as  described  above  oust  be  very 
complicated  and  rare,  since  I can  not  find  one  single  reference  to  thee.  I was 
able  to  derive  the  particular  subcase  which  results  vhen  the  two  dashpots  are 
described  by  one  single  viscous  constant  n and  both  springs  by  the  same  two 
elastic  constants  X and  u • The  resulting  field  equations  ir.  this  still  very 


5 C.  L.  Navier,  (l82?)  Memoire  sur  les 
corps  solides  elastiques.  Mem.  Acad. 


lois  de  I'equilibre  et  du  convenient  deo 
Sci.  Paris  7. 
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general  case  are, 


iv)  ».  [l  - «tlt2]  ,2“  + [\  * ».  - t \t)  ?(?  • + 

+ tlH  {(X«+!  •“>  - »“}+  we  C2  i{’2"+Tf(’  ' “>}  ' 


where 


is  the  Kronecker  delta  equal  to  one  (zero)  for  t^  - (or  t^  i 


Further  t (or  t0)  equal  n/u  . For  t 3 0 and  t0  3 n/u  (here  n 3 U ) , eqs. 

x<—  0 1 V 

(iv)  reduce  to  the  Kelvin  model  equations  (i).  When  t.  = t.  = n/u  (which  amounts 
to  setting  28  3 1/n  and  -3a  = 1/n),  equations  (iv)  reduce  to  the  Saxvell  model 
equations  (Li)  with  3a  + 28  3 0.  The  quantities  t1  and  t„  are  the  retardation  or 
relaxation  times  of  the  Kelvin  and  Maxwell  models  respectively.  The  constitutive 
(ie,  stress-strain)  relations  of  this  viscoelastic  model  are  also  given  in  ref  (6). 
I know  of  no  viscoelastic  boundary-value-problem  that  has  ever  been  analytically 
solved  using  this  model,  which  many  agree  is  more  realistic  than  (i)  or  (ii), 
since  it  can  describe  wider  material  behaviors  and  a wider  variety  of  materials. 
Since  the  "standard  viscoelastic  model"  is  so  hard  to  handle,  we  must  realistically 
conclude  that  all  analytical  viscoelastic  problems  we  are  bound  to  see  solved  in 
the  near  future  will  be  based  on  either  the  Kelvin  or  the  Maxwell  models  of  field 
equations  (i)  or  (ii)  respectively. 

Occasionally  we  find  a reference  in  the  literature  which  contains  a very 
complicated  network  such  as  a dozen  "Maxwell  elements"  in  parallel.  Such  a com- 
plicated netvork  immediately  implies  that  this  is  not  a continuum  field-theory 
approach,  but  rather  that  the  body  has  been  replaced  by  twelve  coupled  domped- 
oscillators.  Work  of  this  nature  happens  to  be  mostly  experimental,  chemical,  and 
containing  little  mathematical  analysis.  These  complicated  networks  are  basically 
intended  as  pictorial  descriptionu  without  much  physico-mathematical  discussion 
of  the  response  of  the  twelve  coupled-oscillators,  which,  per  se,  is  far  from 
being  a trivial  problem.  Sometimes  one  comes  across  ar.  entire  textbook  dealing 
with  various  aspects  of  viscoelasticity  without  a single  reference  to  the  contin- 
uum approach,  or  the  field-equations  for  the  viscoelastic  models.  This  old  fash- 
ioned tendency  is  out-dated  today.  Mechanics  of  deformable  media  has  become  sore 
highly  aathematised  nov-a-days,  than  ever  before. 
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II.  SOLUTION  OF  THE  FIELD  EQUATIONS.  (Tables  1 & 2) 

We  will  now  present  a way  to  solve  the  field-equations  of  linear  dynamic 
viscoelasticity  (i)  and  (ii).  The  technique  we  will  follow  is  common  in  other 
field-theories  {ie,  electro-dynamics)  but  we  believe  it  is  novel^in  viscoelasticity. 
It  consists  of  introducing  scalar  and  vector  potentials  0 and  \jj  such  that  u = 

V<f>  + x and  ^ ^ =0.  By  splitting  the  displacement  field  u into  irrotation- 

al  and  solenoidal  parts  in  this  fashion,  it  turns  out  that  the  field  equations  of 
each  model  (i)  or  (ii)  are  automatically  satisfied  provided  that  the  scalar  and 
vector  potentials  satisfy  certain  scalar  and  vector  telegraph-type  equations. 

The  vector  potential  has  three  scalar  components,  but  since  the  solenoidal  gauge 
condition  v • i|>  = 0 must  be  satisfied,  only  two  of  the  three  scalar  functions 
are  independent.  Those  two  (i|/  and  x where  p is  not  $ |)  together  with  the  scalar 
potential  P,  form  the  three  independent  scalar  potentials  that  can  be  used  to 
solve  the  problem.  If  these  three  independent  scalar  potentials  satisfy  three 
scalar  telegraph-type  equations  then  it  can  be  shown  that  the  field  equations  are 
automatically  satisfied. 

It  is  then  possible  to  expi'ess  all  the  stress  and  displacement  components  in 
terms  of  these  independent  potentials,  which  are  determined  first  by  solving  the 
telegraph-type  equations  they  must  satisfy.  In  this  fashion  we  can  determine  all 
the  stresses  and  displacements  in  the  body  needed  to  completely  solve  the  problem. 

To  discuss  these  solutions  for  the  Kelvin-Vcigt  or  Maxwell  solids  we  have 
constructed  two  charts.  (Tables  1 and  2).  Table  1 deals  with  these  viscoelastic 
models  for  the  general  case  of  arbitrary  time-dependence.  Table  2 covers  the  im- 
portant case  of  harmonic  time-dependence  of  the  form  exp  (-list).  The  coordinate 
systems  covered  in  those  tables  are  the  cylindrical  and  the  spherical.  The 
cylindrical  system  is  studied  in  general  in  columns  3 and  C for  the  Kelvin  or 
Maxwell  models  respectively.  It  is  also  studied  in  the  (t-independent)  plane- 
strain  subcase,  which  13  applicable  to  infinitely  long  cylinders,  in  columns  D ar.d 
E,  for  the  Kelvin  and  Maxwell  models,  respectively.  The  spherical  system  is  also 
presented  in  full  generality  in  column  F for  the  Kelvin  model . It  is  also  given 
in  column  G for  the  oxi-symaetric  case  without  azimuthal  dependence  $,  again  fcr 
the  Kelvin  model.  The  Maxwell  model  is  not  covered  for  spherical  coordinates  in 
these  tables.  The  basic  result  of  these  tables,  particularly  Table  2,  for  har- 
monic time-dependence,  is  that  if  we  introduce  independent  scalar  potentials 
which  satisfy  the  Kolmholtz's  equations  with  complex  propagation  constants  given 
in  item  (10),  then  the  field  equations  giver,  in  items  (7)  are  automatically  satis- 
fied for  each  of  the  cases  considered.  Furthermore,  the  displacement  and  stress- 
field  components  are  found  from  the  potentials  by  the  relations  in  items  {«)  and 
(8)  respectively.  Notice  that  the  complex  propagation  constants  giver,  in  item 
(10)  are  related  to  the  elastic  and  viscous  constants  of  the  viscoelastic  material 
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4)  Displacement  - independent  ‘’qtential  relations 

*2  ~ V(<4>  r $)  - t,  Va  V ♦ a Vip  ■*  ^7*[(,X?i*0*l, 

OR  N component  pqhm  : 

&lv*  SSI-tbC^i)  . 


MEW!  : 1(1*  • ,»,*)  , 1(1  ••  \p£r, «,  E.t'i , X«  X(r,  8,  »,t)  , tf  - 

note  : IX  * 0<<>  . t* » v ( $ =o  . 


^placement  relations; 


5)  stress  - displacement  ^a.AMONf-  • 


V *SA  * §§  S * w*,««4  A*  u*.*-  * 

€y  - £ [ U * Uj.,]  as<a,0>. 


* 2 “ ut  «i»«rS 

ARE  GWEN  IM  S)  AMD  A - V*  U • . Uf  - V YJ*  * S^‘ 


6)  a 1 HA»W  : NOg PRUDENT  t../  EhT-AL,  ^ECRTiONS 


**>  .IRAN  IN.'f  l«E»(!*  '» ' Polj.N.tAl.  .•iRuAi'ONS 


HO  WSMSRAC  SAYRStSiOS  CAN  SESJVWM  VAUO  U*  AMY  ZOOVX  SYSTEM. 


V * l^iilv  * - §::(  F Ss(^‘ ^)]  * * ir&  ■ ^ 5^» 

«-• - 

^ * iWl ' Ms ^ v*v)  * ^ JWfe 

i«T*  - a Hi  - 

it*'!  ii  * ( & •■  >’  h - k " * I 


V ?ijL  Tl 


v -TOR  The  OiSPtACSMEwr  Eittfe) 


tak  »25CP*0)J  0<$-Ei)  • a(3k*26)  aS(3R  2S)j>a 


a C2<a«x  * S if*  io2]^(^*&)  * ao»M*s  s)$«($.  C)-asfl«t*2s)j>4i 


v , . sV.  ! 

i>:  [}!<,•  f*«  VaU  a (x**k»  **J  ^ 

OR,  >M  vOMTOHSMT  t‘OM M ' » T v\  ->V 

*tiA*/*»  H * f|?r 

WMSAfi  A *ta  V9  *S£  6IV6A}  ajvovE  IM  .*)  Anti 


A«?T«t .»  a) 

a,  R.S,  P ARC  OAS RATON S QtHNtMH*  OMT*  MOORC. 


«)’.  [ A.  A 2/U  - (A, *2^)  (TJ.ul -(^‘K-^CO.ul-f  0 


These  two  forms  ;>  ahp  are  eauivAiRMT.  v 4>vsm »n  (*>. 


-iB2LE£liE3sfll^.  ES3JLBTj*k  H.ES 


WO  OCMRRAL  E1ARESS.ON  CAM  tS  QYCM  VAUO  <M  AMY  OOOR&  >TST£M. 


8)  fcTRe&'i  » ••M.CKPtxrNMT  M-.Ai.  RtlRTIOMS  : 

5^)  - Hr(vSa^S4] 

wtvMsfrV * i&vSpCvfti^' 4] 

t*»- [A,  • 1*  A 2yM  »/Ar  ‘ *)J 

C^- l/Ae  ‘/W fell  2 • IS)  - fer (V9V)  - ir  £»<VX)l 

tA8-W*/*»fe]l2  rif  ' V**  1 

these  are  Eat  (s')  wux  a»v*^>  am&  s,  as  given  :s£ava 


rrn^ 


APU  TY»> 


WWRN-tm  E?iu  Typ,  r«»K  i^cAjLjvtCTCk  ivitcntia 


(A£S.<l.*3STP»2Q3lV,ffi  *»  2S  C3R«  2s)  P i!iSi 


!.  fl  am).I  V*tfl  -t  iv  • M* 


L)  _ JORP'V'”'  ■'  : -4W  ■ Nt 


1>)  ft  i.VIN  V IV 


1)  DEFINITION  ot  IMF  h^anE-  - r R A 1 n -„a  p.  : 

THIS  SITUATION. WHICH  lb  AFFKOt’RtATE  T0k  NFlNlTELYl 
THUS,  HERE  THERE  IS  NO  A • DE PEN OENC E AND  U,«-  UrLr.9, 
HAVIN0  <p  - <£;>■,«,«)  , Vr*0,  -O  4 * V,  C 

HAVING  if » (f  (r.  6,tT , ip-U  t,  X-X(r,e,t  ' NOTE.  Tl 
L f*E  FORMULATION  IS  INDEF'ENDENT  OF  t+l . ANT  OF  THESE  I 


i 

fr*‘‘ 


30  CONS  T T U rive  St'LAT'JsS  . 

here,  (?  = «,  S I*  * Q-*/Ae|e 

fV 

Clearly,  K -/<«  , /UH-o  v*. *t , ^.o,  A,  -,i«  , vkai 

• A0,a^«0  , t i*Q,  A.O.f  , a,«0  VR*i. 


^ S' 

M.  *•  . r 


•A*  . 4 - V < «.  ■ 


TSB 


2)  : 

; C,r  “ [■ A.  - J , ij  j A * F(f<«  * “»  gj  ) 

I L*9  ~(A«*  A»  - 2.  (f<t  LyUu  ^)f  «« 

1 v-,.  - [A.  . A„  ^}A  "ERE  A-  €,  fe« 

^<9  “ [f*  *u-t  if]  2 Epg 


^ • 2 • * 5a-  * -€>«a  kVt-&-  £ 


[■ 


I 


*1 


PortAITAL  'DELATIONS 


(r7*X)  = V'x]  * ^ V U*[^*(V«Y>  . 

f-  *#• i * &<V*>  . u*“  - **V 

!»•  e\ * 1 ‘ *<A  9,»,t)  , Tf  • l,  gy  * £»  i,  ij  * it  |g  , VJ-  • ?*  SJA  * 5«* 

■ I <t-0  = o . 


S>  STRESS  DISTH  ACEhENT  RECATlQNS  I 

N-^gil^^v^JV  ' V’» i - 

* if  AS£  GIVEN  n 3)  AN0  a " fib  » * 7 ^ 2 


5)  V'?ES;j  - ^ ~j  i’v *'*  KLiAQwt 

1 ^ ^ &]  ^ ^ vn  * ‘w  i,, ; 5i' 


?■**•  i,'i«  » A ‘ • UV  ifX  fA  » Ir  Vj=jf)  ’ *»•*  l\ 


A c •*  * * * A*  £$  J * Uv  ft 

[t^«  ■■  H.  •»  H..  ij  j i - “fN|  . 


•'■  ■ 5*fA * ^ jis;#  \ vix^  * i!a  ^ 

|gvv8.t)N  ^i*  r v>x}  * 

?.  . . ...  , tOUATIOWitt) 

t u ■ k [ -y  W •>  ] • , , ^uJo>- 

♦ Sfc- *•;»#*•»*»}•>&<***) 

U*'**,'4,£,!w,x>' 


*••  •fc(»G8.<Vo|  . 


r<  c<vh6«  * trains  ( «M  , fi«,  , Cw)  VANIAS 

IK  THIS  VAS6  . Th€  ftetjULT  .S, 


*«“  i'%  ‘ ? ' V>  Hi*  * ^**  * • v 


THtSt  RE  CATIONS  AKS  V'o-WNO  »Y  AVOSTiTVYinO  f4}^J 


♦a*)Jt]^«.uT  »f  $$  ;>  [< . Ue|v4  * ^ pj*  HtvsC  . 

.„.(1l,(i..<»  <■*-*) 

^ w«l»r  ' ’HI5CAKTW  TIPfiNKO  UP  *4TO  COMPONENTS  0Y  means  OF 

««*  v*u  - [vV -%'** i - * Q'l ?»  f, 

* 5^  ^ * ]r  J#  * |-§  £«  i A A&  GtviM  AftOve  «N  S'). 

..n*)a*0iQ  ^ 

taviVALENT.  TJ*  4.N5H  IN  («V  THESE  ’’VO  FO«H3  i5  AND  ii>  Aet  VQLKVACENT.  V*  GIVEN  »M  (4V 


f f*A  */*»  |»  K V‘tN  - ’%  ■ T*  ^j*  V ^ 

CN  *^t  & Uv«s ' ^ • h 


wn6«€, 

Tut  Mi  CON*  FORM  <V)  sS  NOT  AS  USEFUL. 


9)  STF^CtS-  'NDfcTn  SpFNr  TOT  E NTT  AC  RCLArjQSS  : 

t«»  ^\-<aV*V  ‘ V fcX** *§&>  * Ir 

A.us4»2^[^(<pa«5). 

<Tn* -«./». ^5i |£) - v^’v) - -V is» 

«*.-  M.^(£ 

^ v gr  N Ap*.  §j,  r v 

TstsE  am  cat,  lst  with  a - v ^ akd  f,,  »s  ss-tN  ;weas.  ,<»>■ 


opcr»^8^w7rrg»,ft..».i— 


8rK S*  ■ &&  is  <* 

• ** Si  1 »|p(- 


( £s S($r)  • ;^A-r  $?  - 7*  .> 


tjt  «av«  “**tr-E  nflaTionS  <s  terms  of  and  ^ ': 
HUTcACE  v*  » 0V  - V.  . v‘  If?  Given  AfcOvC  M 


f“E6G  «ME'  t OS  {S^  WITH  *<i!»  €,.  AS  favtN  <N 


r..r<ATTrt  iTMAfl  1 '.vcrcM  Pi.ircNnAc 


IN  CYLIN-  | L ) MAX  Af£,LL  OElAbTiC  MQD£L  (S£Rl£S)  >N  CYLINDRICAL. 

COORD1NA  f l r7  FOR  Th£  1*1  ANE-  STRAIN  SUBCASE. 


F ) KELV.M- VOIOT  VtSCOEL  ASTlC 


LONG  BODIES  (*,  CYLINDERS'), OCCURS  WHEN  lt*  = 0 AND  gg  CANY  VMUABLtt)  - O. 
fet£  IS  NO  l- DEPENDENCE  AND  It, » TL-  Cr,  6,  f>  , Ia»»  UjCt,  •,■£).  IN  TERMS  Of  Cp , Tir,  , V»  , PLANE -STRAIN  AMOU 
ftr.B.t)  U.o  uv-o  1 4L  - W.cr.e.f).  OR.  N TERMS  OF  THE  INDEPENDENT  POTENTIALS  (f , V,  X,  IT  AMOUNTS  TO 
t(n  0,0,  Wr.s.t)  note  that  sv  eas.  b.ii-b,  V.- -vJ^X.  So  in  plane- strain  only  <p  t appear 

tfiON  'S  independentof  V-  any  of  these  equivalent  definitions  yields  . G,t  = o , €,r  


1)  RELATION  BETWEEN  [HE  TRAC 


2^*.  w t.^3** 


C),l  « BULK  OR  HYDROSTATIC  ST 


l)  CONLTlTUTvH  RELATIONS-. 
$■’’=  Mtt  -zue  trr  } * N, 

^09“  *«-A  + 2-,We  G««  J *■  N' 


M,{AtA}  * N, 

,^TB*  *4  ^r6 

WHERE  A IS  CiV 


C HERE  M,  AND  N,  ARE 

The  same  operators 

SIvGN  IN  C.-8T  OR  E-8)) 


where  a is  given  below  fie,£rr-  ■»£,»)  in(bQ. 


C EvER^  STRAIN - 
COMPONENT 

2 .tr„  , -L  , 3U«  _ Ut  , Of  2.£»h=0  , €*8  ” O WITH  SOC'NDEX 

r 39  «r  r ^ , VANll.nES.  ) 


2)  CONSTITUTIVE  RELATIONS'. 

t)5^'2Cu« 

f^r  = aeAAff|-)AT2^*ME 
1 tH=(X*  + + a C/M+r'v  If 

2^/'A+/V|{)Gr<> 


3)  STRAIN-DISPLACEMENT  relat 


(■'  — 3Ur  , £ = L SjA#  +.  Uj*  , €.-  a 

nr  ’ c»  |-5A  / ♦♦ 


1^1"  [ - ♦ ^-L^iVx  ( i.L)\]  OR  IN  COMPONENT- FORM, 

Lev**,)  • IN  TERMS  OP  <£  AND  X,.  or, 


4)  DISPLACEMENT-  INDEPENDENT  f 


. T-XCoO.O.  v - jifte,?!; 


in  TERMS  OF  Cf  AND  Ip,.  ( ip,  “ * ^ 'C  ) 
v2,=  grY-T-l?  ‘ 


5-J  STUt  Sa  - C SPLACCMENT  RELATIONS 

M-UA  - •iNC^r]]  - TO, 
te«  = M,  { Xe  A *-  ?.  H*  [ ±c  . L j * N, 
tii  - M,T  > t A ( f N, 


EBasareHsaBagBa^^g 


(HERE  M,  t N,  ARE  TTC 

Same  operatoilS 
GIVEN  IN  C-8)«rE 


«:tl  »(*.***&)  A-ai^A/A 


CRUMP  ENT  POTEST. AL.  KL-ATU 


( THKSS  AWE  bOScAT  SOBbTlTotEO  INTO  EAS.C3).) 


WHERE  the  e^  are  given  above 


■ -r  isC'7'1^)]  , £«*  = -p*  fl§k  J + . 2«r#-  2-^e(-^)*C|^--VfF-r'j»«s)  ylx 


STRAINS  ( e„  , €,r  , VANISN  . TO  have  THESE  RELATIONS  in  TERMS  of  AMD  V.  , we  REPLACE  V*X  BY- 

Tse  RESULT  IS, 


tNL.  -o  '-/L-  I I W ^ 


TIONS  ARE  FOUND  BY  SUBSTITUTING  EOS . 4D  INTO  EOS.  3). 


w, mi 


7) field  equations : 
i)~ 


T)  FIELD  EQUATIONS  . 


|:vs;l  +[«+ 

pTORM  VO  is  NOT  AS  USEFUL.  $LV“  ^ r~ 


WHERE  A AND  Va  ARE  GIVEN  ABOVE  IN  < A'l  l CS).  fo*  -/<*/?•) 
FORM  AO  >&  NOT  AS  USEFUL. 


O t/*«VN&JV*4  ♦ 

' [/*«  */*«  - FTwT* 

_ t/M*/w  JJ[ v*t4  - >^f»r  ♦ fIs* 

it)  [ x,  ♦ a^t*  ♦ ] \)( 

THE  DILATATION  A IS  OiVKN  IN  * 


1 

4 

"HR 


E^EMPg-lNT  P^OTfcNrtAt.  ^9- 1 AT  I O N ^ 


4i*i 


r||]  V’rp  t - $?(?  Is  ^V*%))I 

|f]  Vacp  -h  2 [p, + r- lr  + lp( r Is ( v”^)] 

rfdV*<P 


6)  STRESS  - inDEPP.NDEKT  POTENTIAL  RULATlONfc 

f«rf  » Ml  i !V  * [i£S  - |?(f  CY**|]  ] ♦ Ml 

1 aH-M,lW7a(p)*  N,  . t«,-0  . Cw-O. 


8)  STRESS-  -NPEFENDfiNT 


iBEtkttB 


SE  RELATIONS  IN  TERMS  OF  vp  AND  V*  WE 
X BY  - ty»  • V1  IS  GIVEN  ABOVE  IN  C?). 

S.  CB)  WITH  Ai  » V»CA>  AND  €(|  AS  QIVEN  IN  EOS.  (.6). 


(TO  HAVE  THESE  IN  terms  OF  ¥ E V«. 

* we  REPLACE  ‘Va'X  BY  - 4>»  ) 


r jt  ^ 

N . feLl  7 > 

^ R^^fTjL  ♦Prr8?  J 


i«-“  CX**  Xr  |fl  Vss)  ♦ -fv 

Cf*  [/*.  *•/**  If  Kv&»  'V  fc(* 

cv*-r^./vyt£lr\^SaC<P^ 


' H -ToTJi  JJ 


,1Y  ~L',Ev '<A,  ► LlY  I UN  THE  SCALAR  A VECTOR  PQTE.N1  i»LS 

( 9 ft  » 4w  Si  t DO  oVw  • k 4I4BF3UF L\  & ♦ 2p(a« *aP)  ll 


9)  TfU  ...RAf>H  F.QUA.'JN>S  L‘IR 


O 0*rt|r)y*«-4l$R 


frlSCOEXASTIC  MO  PEL  C PARACi-EL)  iN  q£NER,AL  SPhCRi^AL  l.  <JQR.D  I NATE^ 


EEN  I HE  TRACES  i'h  AND  Lal  Of7  T~H£  ^ ?Kfc.U>  AND  STRAIN  TENSQRja  ■ 

; ^ - [C3A..  * 2/W)  T (3Ar*  2-fr)^]  tfcfc. 


0 t^E'IN-1 1 


HYDROSTATIC  STRESS 


E RE.  LA  T ION S - 


*.  cubical  dilatation  »A 


WH6R£  A”£*<E‘  err  + W£*A  “ U*Lfc  °R* 

• . _»  riv„  / >-  1 A M.  2 C LL.  *.  Lr  ^ 


ijflk  T 2 e.M«  + L rt. 

■A  + 2 fyw*  Mir  ^j)  fc«o 

i>«* 


*5,- 

“ 2(^4  ■‘•^ir 

tr«- 


2")  CONSTITUTIVE  KE^ATlONS  •■ 

V - L^*  ^ y A + V**&)  *rr 

c«« - [A  + ^Ie]^  •*• 

?>4>  « [ A*+ A + 2 f|U*  ♦ |e) €*♦  » 

Ltr8-[)J«4^u|-]2€r9  4'V»VSi 


Aim-  DISPLACEMENT' 

RELATIONS:  ("2E*.  ■■*  Efe,* ‘O') 

■ €«♦  - Ifr  «■  ¥•  “te  » 26^=  T jSr-  ff  ^ 

0 


| <^A  ♦ ♦ A»^-a-j .,  .- 

jlUE  GIVCH  ABOVE  IN  3)  [N  TEfcMS  OP  THE  DISPLACEMENT- F1E1.D  COMPONENTS. 


NDENV  OOTENTIAU  RELATIONS  *. 

♦.gF-l>VaV]  = TllTr^tr^'^lS-  • 

iS>T§Tr>'W\|  ■ VJV  - V ^ Si*  ‘T*1^  *'T  frikp  ' '^1  srS  * 3^ 

-^T5  5555'. *■  • 

♦ 7rfe*e  -^S^S  ' ^ ■ 

$ C*i»“&wKva*) 

»■,«  ,.  ....  ..„  t v.  .«  «»0  ««»•  AS).  


r*d  ♦ W&>  -P  w °*  lN  component-  porn, 

8^  - V#  - ^ $Sff]  - l*.  *P* -<A**N&  1 # - P 

^[Vi*^*CV*MgrJ4rJ|  -p  §>  m 

$H  *IA  •/*  ^ * i 

«•>>&]  Wa>  •f‘**Md**<*'a’  - 

l»  OivtN  IN  #)  ABOVE  AMD,  V*«  & ♦ ^ |p  • E & • $4  - 7*&*$  $p 


CEM&NT  RELATIONS 


SUr 

Lf  $f ' * VJ  r ^ as  im  ▼)] 

[£.  ^cotO]  L J 


o')  ?T«AiN-  Nn£PKE4DCRT  pqt;\TiAl  RLuATiQN S : 

eA*'l?^r*K ' VX  V - >'  Sr  * ' 7 ^»)  ^ ' 

2€m-°  , 2fe^-0, 

2.  |?[-  ^ (lip*  gp.CTV5)]  - U)-  ^ ^ - %» S8-  * 


?)  -\S.lQ  giauA-'Ona  -. 

*•)  (/‘A*/**  jyi<u*-s$f  ‘fii'if- 

#fcwl 

WMfitt,  Vj  I, 

R3«M  <»5  WILL  HOT  te  6IVKN  MtMJE. 


P6HT  POTI-NT'  *L  . -■vE'l.C.  VlQN 


9)  bTKCOS- ‘HCftfE-NOSHT  IVMT.NTlAt,  R.ELA1  'ON«  : 

«Vr.(A,.A»gt)VSl4  2(f,  »/Vgt)[|*-%(«  * g?<T^)  “ ^?Cr  W«V)J 

WC*«*  <*•] 

<A>a-2(^.»P* g«XC+  If  ' &)5 V • S? c*-^ - vi v ! 

s'°  ■ er,'°  ' 

I.IFiltv*  fll«trV)))  - ■ V)tv‘  tt]  (nmWM  ALn.OyO-  yT..^  9U  Ls  Nor  C0Nv6N«*r  .O  I..M.E9S  .4.nui* 

«S»*5Tin«rTO  TH656  «KLAn6NS!WOH  s*4)  .H  TtRMS  OF  ^ AMO  v 

^ (V* v) * SW9 3l(ir- r)Cv* a:)]  ,HTO  c«*  a^-i  (-un-g  aac  t«8.  •»  •.•«»!. nrurso  jhto  eo6.  aO  (A-  5^o>). 


y*p  * ity* */•*&}$»(«>  ♦1?'.'*^)-  j?(r  v'v)] 
v+zip  Y*  yKVl?  * U &X*  * * v*v  - -^  y^»  k (w,;°)] 

- 2 [f*  *r*  hl$k*  w b ’wSw*  ^frV))‘  ^‘v  4 &wi»  l»  • 


■ ,;  ULhiiAPH  LOS).  '-  CK  ]LA1  fiK  A VI'.L  fOM  EOTT  NTIAlA 


i)  (i  ♦Hjj)  - j||  ^ . i»>  0*N 9)  • 


at,.-  £ j£g [♦  * $8] - 7 §s < W)  + hr^^  “ $ jS& + rfeti- . 

2 ^elrl^  + 1*  11  +[a?5“7lr"'Kle,‘l<'Va^  “ 2 tr [?!!]*■■ 


2.6, 


7)  THE  FIELD  EQUATIONS  : i.  FOR  THE  DISPLACEMENT  FIELD) 

»)  '■ 

£qr+2S(P*Q)]9(^‘U)  ♦ <2(3R+26)  ?*£■»  2Z(3fi+2S)f>  u. 


uy. 


2 i[2aR,  * S (Pf  2Q)j  v(^.  U)  - a C5RT2  s)  C)  * 2S(3R+2S)ji5 


note  : v*ii  - v(&‘U)  - Vx  (V*  &) 

a,  R,S,  P ARE  OPERATORS  DCPSNOINfl  ON  THE  MODEL. 


7)  FIELD  EQUATIONS  '. 

i)  ■ •*•  ■*•  C^v-  +^U.^  1 ^ %£T 

or,  m component  form: 

Aklly4u«-  Vt> 


*[*.♦/«.  ♦cw^§r]f|  - f & 


where  A and  Vs  are  given  above  in  a")  and  e), 

~v.  [ A*  y 2/A*  4-  a,v2^)  l^lv^.tLHpe^lg^xC^-f  0 


THESE  TWO  FORMS  i)  AND  U)  ARE  EQUIVALENT.  V GIVEN  IN  (4\ 


7;  field 


o IH 


TH1S  CA" 

V4u  - [va 


<*■'/  [l  • 


♦5*1^ 


THESE  Ti 


0)  S.I&ESS  - INDEPENDENT  POTENTIAL  RELATIONS: 


NO  GENERAL  EXPRESSION  CAN  BE  GIVEN  VAUO  <N  ANY  COORD.  SYSTEM. 


&)  STRESS- INDE  PENDENT  POTENTIAL  RELATIONS  .' 

^Vp*5  [2*Y2v|£]ya^*2(^*^|g)[^1((()  + if)  ~ |?( 

%•-  ♦2<MV  &>$»  £♦  & :w»a3j 

*»-[*•*  2. 

*w- ^*^Ir  - Wv*  V)  ^rcvx)] 


S')  STRE< 

tfr,.-  pU 

Ci»-M,l 


£*<**%)  - - V £»(*•*■>] 

tA*-*[^+/v|{][2  |^Cf^+^i)+rin-7lf  *^^Ovax] 

THESE  ARE  EQS.  (£)  WITH  0 « AWO  €.,  AS  GIVEN  INEQSfib' 

CaB-M, 

" 4F*- 

THSSE  ai 

<J)  TELEGRAPH- TYPE  EQG.  FOR  THE  SCALAR,  t VECTOR.  POTENTIALS’. 

0)  TELE  6 

0 tap 

CQ.  ii) 

t7! 

*)  [4Q.R.  ♦2Slf*2cO]v4<p  - 2S(3R*  2s)  ^ 

- CL  [ V * ($«  V)]  * 2S  j> 

Ed.  ii)  CAN  ALSO  BE  WRITTEN  AS  FOLLOWS  : 

v*  0 «■  as£  and  'tf.ij  «o  . 

*> 

it)  *>[l*  “ qJ  * N*' m ‘ 

EO.  ti)  CAN  also  be  written  AS  FOLLOW*: 

0Nk3  V*v  b ^-V-O. 

■M8Mg»EBra^^ 

1°)  TM 

i) :  [2|S* 

ii)  V7\ 

FOR  6}' 

V*i 

a 14<2r*23(.pa2ca]v,<9  - 2s  (3R*  as)  j>  ^ 
<«  a v°v  - as  $»  ^ 
iu ) a v*x.  • ?.s  $»  ^ 

i)  ^ ^ , 

it)  [1  * N Jg]  D4v  - £»  . tJi)  [l  t N If]  - Ij 

' N’^  * °*‘*-*f^  > 0''“f  • 

«)  SQLENOIDAC  SOLUTION  OR  VECTOR  T6  lE  GRAPH  - EON. 

II)  SOLE NOI DAL  SDLUtiON  OF  FHi  .KC'OR  TELEGRAPH  - EQ 

UATION 

THE  SOLUTION  EXISTS  ONLY  In  FIVE  COORDINATE 
SYSTEMS  (.CYLINDRICAL^  SPHER(CAU)  INCLUDED) 
BUT  IT  VARIES  with  ThE  SYSTEM  USED,  AMD  NO 
EXPRESSION  EXISTS  VALID  FOR  ALL  SYSTEMS. 

NO  GENERAL  EXPRESSIONS  FOR  RELATIONS  t*),(4)ORv6> 
above  can  BE  given  for  this  reason." 

C«EF  I P MOM*  S W.FtSHBACH  4 METhOOS  OP  THEORETICAL 
Physios"  VOL.  2,  CHAPTER  A*.  RcgAaw.hilL^ESS) 

It  CAN  be  SHOWN  that  :iiE  SOLGNOiCAL  (>«,  9-  o •o')  SOLUTION  < 

ifiS  ) ist«e  solution  of  [Un^-]  ’vi*v^»ri)  + » *»°  amo 

tJ)  - Oa  Cv  t»)  * *9  x C^«  <6»  *)  • C^s))»  t,  - ^(||) 

WHSRK  VCT.e.A.f)  AND  fc(r.  e.A.t)  AKE  Two  new  scalar  functions 
60UATIONE.  W«.  [t  . N |g]  V14  V - LJ  AND  [t»N|g]Vs2»^ 

• 'r  * '•V*  ■ $$  * V * vV  * - 

WHICH  RSLAVe  THfi  THREE  CYLINDRICAL  CpMPOMCNTfc  Vr.  V,  . V»  OF  T 
V AND  X.  SUBSTITUTING  $ FROM  THE  SOLUTION  (.A)  INT  O £ f» 

cylindrical  components  of  the  Displacement  field  and  the  sc 

AS  UWCN  ABOVE  IN  EQS.  t*).  “**  SOLUTION  (A)  0«(J»)  .NV-OLVES  ON 
WAVi-EO.  ( OR  TELEGRAPH. EQ.  WilHflUT  ITAMPtNO)  PROVIDED  THAT  1*1  A 

tale  the  compartments  To  the  Right. 

HI1MAUK:’.  : 

a.)  The  BASIC  RESULT  OF  This  TABlC  , *N0  That  IS  WHY  IT  WAS  ORGANH6I 

sati&figo  in  each  cine  or  rue  cases  considered  and  then  the  tmi 
BOUNDARY  value  PROBLEMS,  and  which  S common  M OTHER  DlSCII 

re  HEADER  TAMILIAR  ViTh  F'^CtROOY  VAf.c  T|,^oRy  WILL  QUICKLY  N 
b)  it  IS  SVIDENT  That  in  ALL  thE  CASES  in  THIS  CHART  WE  RECOVER 
C)  IN  *.VLIVDHlCAL  COORDS.  TuC  PcAHt-  STRESS  RESULTS  CAN  RE  QUICK! 
X ) THE  CASE  OF  Plane  - STRAIN  {OR  PLANE- STRESS  IN  VIEW  OF  C)  ABOVE)  Wit 
© A CASE  OF  SOME  INTEREST,  WHICH  REDUCES  ThE  TWO  VI6C0U3  CONSTAN 

!ft"M^-V,'l’}  4 &<***! 


MTOCQ&CJO. 


ar  • 


* v 


..j» 


7)  RELP  EQUATIONS  : 


equations; 


i)  L1  * a ?tlW ’*&  - - 5»  ^ ] + If  [vJ<i  * 

+ 4r  v(^.  U)  - 2» £$]  - gj  H • (°»: 


THIS  CAN  BE  OPENED  UP  INTO  COMPONENTS  BY  MEANS  OF 

Sr  * £ H £}  ♦ & «•  * A AS  GIVEN  ABOVE  IN  «). 


[j«»  ‘/V^][vsur-  %,-  ^ VV)Jt j ‘i*  ?& 


* U (v-  u)-9*  iV* u)  - ]-  5*  [I?  + ^5*' 


WHERE,  A-  *!,+  <!*.+ 1-^  , V-1^4!-*^ 

THE  SECOND  FORM  ii)  IS  NOT  AS  USEFUL. 


f),  TMF.SE  T^O  FORMS  i)  AND  ii)  ARE  EQUIVA CENT.  V*  GIVEN  IN  (4). 


Q~)  STRESS- INDEPENDENT  POTENTIAL  RELATION^  * 


' - INDEPENDENT  POTENTIAL  DELATIONS.’ 


S')  STRES 


m?: 


frr'  *t{**v\  * ” *< 

(r4\ 

C„-M,tA*v4<1>  + 2^4[^!l((p>-^')-  |5(v1i)))]}«'  N, 

lSsC^^)4fcwH))^fv**] 

- M,/u[a.  - §r  Cv*V)- 31 

<*.-  m^o  ;wiR£> 

M .,  ■ 1*0  ]wV 

^ S*  * (3o*aiS)*-  gjvi.  4pl-  i 

:|HEQS.&|  P^tsE  ARE  606.  <S)  WITH  C = Va<r'  AMt-  C<,  *S  GIVEN  1NEQS.CS). 


'Kr-  V*«f  ■*■  2[^l«+/<v^][^|  - SrC?  &S  Cv^Z))1 

'Zef  [A,. * It ] vaq>  * 2 [/*•  +-/V 4 Ir  + ^4^v^)] 


TO  HAVE  THESE  RELATIONS  'N  TERMS  OF  AND  V«  v* 
REPLACE  v2  X BY  - W*  • Vs  IS  GIVEN  ABOVE  IN  (?>. 
THESE  ARE  EOS.  C5)  WITH  A»  Va(f  AnO  6<j  AS  GIVEN  IN  EOS.  C6). 


0)  ~ELSGRA1N  EQ 


asssHB ess 


BiAUJUHB 


W 


4>  C^^'4i*felVa<P^[(3«vAp  &^*apCK*3pH] 

ii)  ->  (^*0)  « It) 

KO.  ii)  CAN  ALSO  6E  WRi  TTCN  AS  FOLLOWS, 

?*$.•  *ND  V-V*0 


^ • U>OvN^V^  -1  • 

WMERe’  m.  lr* |jjy  , n«  , cf»  , -■/  - ^ , ^-\5-o 

HERE  , V - - -X  V2  XCr./.t)  . HENCE  , 1 <?•  t?  »0  , 

o •-#-  * lb  I -r-TCAl  ik>  Cl'niCC  ii  TfDKJ  ATiwC  FORM 


HERE  , U)  - g.  . -i.v;  X(r.ff.  t)  . HENCE  , V- V » 

AND  £«,.  ii)  CAN  i>t  WRITTEN  IN  E ITRER.  ALTERNATIVE  FORM 
'SlVEN  IN  'O’). 


:SKi 


«.)  (*?♦&)  • M)  v“^"gtaP* »')!?'  ‘ 

FOR  Cj«C  1 (A,  8V68CR')  EOC-)  REDUCES  TO  { 

^1<P  " ♦&«•♦**)  t 'J  H8>  * ei’¥‘ 


*>  C1  * M it) vi^  " I?  * ,HS  e<a 

AO  “kf^v»  0*  , -V^(UN|t>V-^]-»0  -Vi 

WHERE  ►.,  N,  Cf  , Cl  ARJK  GIVEN  ABOVE  IN  3). 


. F QWAT I ON  iN  CYLINDRILAl  <•  .QOWO'N  ATS  & ■ 


•O')  SOLUTION  OF  ThE  VECTOR  tGl6^KAOH  GOoATiON  (<'1,  {<* N^)7’,^< 

- o ano  that  it  is  given  by  e ither  one  of  the  following, 
it  - ^vH')  ' *»  V**’  <*■>  t CYLINDRICAL  coords  only]  , 


ii)  ',.)i  C.nqiPAl  ^qlv/T'Cis  qc  »».(,  l y qr^  LgM^Stl^y&llS&LlLliJlltL 


T*4?  SOCUTtOM  FOLLOWS  fMA’T’  (jtVILN  TM^  t.EPT  FOK  GGN6.KAL.  Cvu^NDi 
P^NCG  » TMGHUTORC*  om  CAW  ANALn^OO^LV  &WOW  THAT  t'M £ $OV.K|t< 


0 - <5  • (Ox  c^Xcr.F.OI]  - - it  v4*- 


tAR  FUNCTIONS  Cl*.’  V ts  KDT  t^l)  SATISFYING  SCALAR  TKlEGRAPH- 
tN^fcl  V*X  m 4 a-4  * the  SOLUTION  Ca)  in  COMPONENT  FORM  iS, 

h-  . . H * IT.  l >** 


#3 

7 ' Si 


r » •A*  • <5 

V V V»  OF  THE  VECTOR  POTENTIAL  y TO  THE  TWO  SCALAR  FUNCTIONS 
i>  INTO  U » <9<v  *•  "&«  55  yiSLOS  ni6  RELATIONS  BETWEEN  THE 
D AND  Vh6  SCALAR  FUNCTORS  (i*  . iNDE PE NO€NT  POTENTIALS')  <P,  V*  * 
) involves  only  6 pace  operators,  scxta^o  hold* for  1 »e  vector  . 

IDGD  THAT  1*1  AHJ  X TIRjN  SATISFY  SCALAR  WAVE  GO G. 


PPO..D&D  That  tl>,«,t)  SATISFIES  Tl(6  S'CiLAR.  TELEGRAPH  EduAh- 
SINCE  The  LAPLACIAN  AND  ThC  CURL  OPERA  TOK.S  ACTON  A a-lNOSPENDtl 
CASE.  , 

SOLUTION  (A)  iN  COMPONENT"  FORM  IS  • ty.  •»  0 , 1^-0  , W*  * - v. 

A *•  COMPONENT  WHICH  IS  »-  INDEPENDENT.  8j6S'"'TUTINi  SLN.  CA)  IN' 
POTENTIAL  RELATIONS  &ivEN  ;N  EOS.  C*'1  ABOVE.  WE  REPEAT  THA’ 


TmE  TELF  GRAPH' COUATi ON  REDUCES  TO  ThC  WAVE'  EQUATION  , AS  IT OCCL 
SCALAR  wave -EQUATION  [V»  - i.|^]  X -O  . 


I WAS  OR  SANDED  ..«  SHEWN -»OV£.  -S  THAT  iF  WE  INTROCUCS  INDEPENDENT  AChLAR  POTENTIALS  WHICH  SATISFY  Tug  TELEGRAPH  - TYPE  COS.  10' 
ND  TM^ThE  DISPLACEMENT  AND  STRESS  FIELD  COMPONENTS  ARE  GIVEN  N TERMS  OF  THE  ROIT-nTIALS  BY  EOS.  A) 

(MOTHER  DISCIPLINES  fit,  ClECTKOOvNAMiCS),  IS  APPARENTLY  iC/EL  in  VISCOELASTICITY,  ano  wife  IS  WHY  This  CHAi.T  was  DE  L 

W.  Quickly  notice  t«e  analogy  of  me  PRepEnT  atio-n  >n  this  Chart  wiTh  that  op  elucthojiagaatic  theory  in  a coulomb 
t we  RECOVER  THE  SIMPLER  RESULTS  of  dynamic  elasticity  in  the  absence  of  viscosity.  t‘*.  Vr-0  .H.-0  tor  the  KELVIN-VOiGT  MODEL  ( 
* CAN  BE  QUICKLY  OBtA.MEO  FROM  T«E  PLANE-  STRAIN  RESULTS  IN  CASES  D)  I E)  BY  THE  STANDARD  THICK  OF  RE  TV.  AC*  NG  L.  IN  PUj^E  **TRA^ 
fiDF  OA«JV6)  WITH  AKIAL  SYMMETRY  ABOUT  ThE  K-AAiS,  CAN  BE  OBTAINED  FROM  O X Q'i  BY  &E  TYING  U,  .O  S ^ [ ANY  VARIABLE]  -O  . THEN,  ThEKE  IS  ONLY 
P<2^0N^ANT6  T^WG.  OCCURS  l£eN  THE  BULK  VISCOSITY  iS  KERO  OR  NEARLY  SO.  THIS  AMOUNTS  To'SaVinG  U..  Zf. -O  (KELVIN- VOIGT)  0 


7)F<£LD  EQUATIONS/ 

-q[C3«+4^+|p)^f  + ipctp+**)  |£]  . 


7)  FIELD  EQUATIONS  : 


[3a  + 2p*-^][VaU«-%» 


WHERE  A AND  Va  ARE  GIVEN  ABO/6  IN  (4"}  ( C.5).  Cc/  m }***/$•) 
FORM  <i)  IS  NOT  AS  USEFUL. 


i5  mu 

' [fc^rfeU  V>Ur  - *Ss«  */- 

tyw*/*. fcltv’-tL* - 71 ^r*7wh^V¥  * £*!&& 5^*i */**  ^ 

<*)  [A,F2/«.»C^-a/v)^3  v(v.«.)  - f® 

THE  DILATATION  A is  GIVEN  iN  5)  ABOVE  AND,  V » |p, , ♦■  2,  §yi 


Vax))I 


8)  ^TgFSE-  INDEPENDENT  POTENTIAL  DELATIONS: 


v1* ♦z[^e»-/v|tK|^(<p  + S?CrV'))-lf'(rVl!*>2 
vlvvklrt^t  kW  *i*63(**  8,CrVSI-v*v-+ 
^.[J4  + Vy  ^2[^^|fKr^"^i5  + r^lpX<<,+  lrCr 

Vt^felWbfrttkfa*-  gpCrv))}  - iA»kc*N0 + Ss($?  -4- 

C/»**/v  y [2  M-k Is C<p  - (••<«)}  - la  Cv“**> '■'hihib-W 


<)  ‘Iclv^p  *^U4P*--5«*y  5^  + 2ff3c‘-2P>  If} 

-•)  A“°  *,i3"0  * 

• for  (ie,  u««Xt  w >)«  Va. , i*,Ruooe(0  c^.c*-)  'S  •• 

i)  v*<t 


..ala*  AND  VECTOR  POTENTIAL 


I)  (l»M|r)Va<p  . iO  (l.N|f)^(«^)  + 

’ Hmfc  * * C*1*/ 

THEG6  EOS.  ARE  OBTAINED  BY  SUBSTITUTING  tt  - 9(p  ♦ INI 

OF  the  FIELD  Eas.  r).  EQc  <U)  CAN  ALSO  BE  written  as  f 

I o'  Cl  • X &)v-  E-jjVft  ^ g'$~° 


TEnTiAu'j 


tenTialS.: 


THE  ea.  FOR  <9  IS  TM6  SAMS  AS  A6QV6  IN  <»).  THE  EQ.FOR  X.  IS  i 


.v«{vU.-^(ap-&)^-o  0,1  'ntsrmsofv*, 

v*v. 


where,  m«  ^rl^i  • N* ^4  » ^ * c*“$ 


agBBBS8BB«gBB8Mffl 


.0O«C/ihAT6'» 


P FOR  GENERAL  CYLiNbRiwAuLOO^CMNATEt,  SLIT  NOW  WITHOUT  tMF  A PEPfe'N- 
HOW  THAT  THE  SOLE  lOltXAL  SOLUTION  OF  (\  + Hjg')  V*  $ “ £|  '&» 

N]  *•  - e,  v1  x tiO , 4 

t,  telegraph  COoAfriON.  \l.».  v'uN^)V*Z-it^i  • *^J**TY_ 

r On  a e- independent  function  , they  arc  also*  a-tndupcncsnt  in  this 


it  CAN  fit:  SHOWN  that  The  SOLUTION  Of  EES.  4)  L*  )*  ABOVE  IS  <h>j 
PROVIDED  That  the  TVO  NEW  SOLAR  FUNCTIONS  Vu-.IMi.O 

..  ._  _ c . <P.  ..  H * ^*iv«  _ \ A?  V Aklt>  Mil 


PROVIDED  THAT  rMt  TWO  NEW  slalhh  -w~v 
soalar  telegraph  eas  G« , 0 * N h ) v V " a at^'  *N0  0*N» 
( NOTE  that  V 13  NOT  lljl).  OPENING  UP  >HE  SOLUTION  t* 


W..O  . U.-  (B)  . so  THE  VECTOR  POTENTIAL  HA6  0NLV 

»TiTUTll*SLN.  (A>  INTO  ii  • * *&»  IJJ  ViCLDS  ThG  DISPLACEMENT” 

WE  REPEAT  THAT  all  THiS  HOLDS  >N  TNE  ABSENCE  OF  DAMPING  WMBN 
Equation  , a»  it  occurs  n Elasticity,  provided  that  % then  i^nsFifi*  a 


(note  that  V 13  not  |$|).  opening  up>he  solution  tA 

v,“* ^ u * ***  - *4w ^ sla ^ * v ft  * %£ 

WMICW  RELATE  The  J SPHERICAL  COMPONENTS  OF  T«G  VECTOR  POU 

solenoioal  by  the  Gauge  condition  to  the  a scal 

SClOTiON  iA)  >5  valid  IN  SPhERiCALS  ONLY  , AhQ  WHEN  1T16  0UDC1 
TuE  RESULT  16  The  DISPLACEMENT  - INDEPENDENT  FOT^ITiAL  REL 
AUXILIARY  RELATIONS  NEEDED  to  PROVE  THIS  are  T«UvU»*i  J * 
^.r?v,xT*  . IN  r«E  absence  OF VISCO &ITY  l« 

CQ  (A)  ABOVE  IS  THE  SOlENOiDAL.  SOLUTION  OF  THE  VECTOR  WA 
PROVIDEO  that  Then  V X SATISFY  i SCALAR  WAVE -SRI 


GRAPH -type  EQ$.  10)  Then  The  FIELD-COS  ?)  ARE  AUTOMATICALLY 
RE5PSCTIVGLY.  TwS  reCKNlOUG,  WHICH  IS  BASIC  TO  SET-UP  AND  50LUE 

s developpeo  for  all  These  casee  and  viscoelastic  models, 
omb  C>».  ^-3>  -O)  rather  than  a loreny  t (.<e,  )?•$  • t * ° > <>*vse . 

tVIN- VOIGT  MODEL,  AND/OR  <A*0,  B-O  FOR  THE  MAXWELL  MODELS. 

THE  PLANE -STRAIN  CASES  BY  X,  DEFINED  TO  BE  . X.  - 2 **/*•  A *•  ' * 

Q.  THEN. ThERE  IS  ONLY  r .DEPENDENCE  (•«.  «,«Rp<p.«  » ljf.(p<P.«)  ••  Uy  * W • 
^»0  (KELVIN-VOI6T)  OR  3w»af>-0  CMAXWELL)  SO, formulas  Simplify  »mT. 


»1)  REMARKS  '• 

«•>  REMARKS  <0,  b),  •)  OF  ThE  CVL.NORlC.AL  CASES  TO  ThC  LEFT  ALSj 


ii)  NOTE  AGAIN  THAT  «N  SPITE  OF  THE  APARENT  GENERALITY  OF  TUI 
ABOVE  WITH  ITS  VECTOR  OPERATORS  , iT  5 ONU  'AUD  IN  5Ph£« 
FURTHER  that  THE  TWO  PARTS  C>«.  0‘l’i-V'  *m0  U«jvALr 

TREY  ARE  USUALLY  PERPENDICULAR-)  CVCN  WHEN  V X . 


i/i)  t»E«K  l#  NO  WAY  TO  HAVE  A PLANE- STRAIN  OR  PLANS- &TRCS* 

Tut  closest  wc  can  oiT  to  Them  is  the  axially  symmetric  a 


■:y.  p;  wh*,^  'TrjffJJJ 


iwn- -?*HS83r  ^7^9S^W  ■ ” ” „r<,  . , , ,,<*  » , , ,*«  » a*  ^ 

i,4*.ft.-**tt  >Wik«)l*W--«»-'b  wm*« ^2l__ 


(jg  ECi.S  (4-)  iM^O  £*<a&  C3}» 


^♦^.a^v)|j.]«fv.tl)  -f>^&  OR  in  COMPONENT-  PORM, 

r ^ £]«■[*»♦/*•  ft  “*&**  ,* 

]vCV-  tL)  -lf*++fl r|t]^*(^U.)  - f @ 

> In  b> above  *nd,  vl- $£,  *rtiro  t?» 


I 7)  ~ie.lD  eauATtONS  : 

<)  [/<**/L,|l][v4iir-^-M3f-  +[^*^(V/^|E]^t“ 

-f&fr 

$»!$■-  ^fcw) + £** ^+a.+)^) y • 

WHERE  , ^-|_l».2:.|_<.-m:5k<-  (j$$-  §* 

A-  |«f  v 


FORM  iO  WILL  NOT  BE  6IVEN  MERE. 


DTENT i At  DELATIONS  : 


*t*kWrk'k&fo*  l?^))-^V- ^WstwifC^il 

*/v  jjlK-Vfc  * ^ & * tto*  * ^)  - ?4v  * h ■ 


S')  STRESS- INDEPENDENT  POTENTIAL  RELATIONS: 

^*(^*^sE^^+2^*f^lt)[^C<i,+l?<:r,*,))  - &*(r y«V)] 

V(3.  *^!e )v  ~ v~  ^ 

**-(A.  *K  fe)  +•  2</Jr  Afr  feXflp  & -■ HP  SlX^ * is" ^)- - v"  V 1 

- O , tr>- O , 


r<p*  ip,crv))}  - wb^(^*v)  * &($f  - V)0* *>]  • c^ese  **E  «*• 

' f « SUBSTITUTED 

k^^W)}-i»^,W--nWI|(ft'-<WCv*x)J.  -to 


C*  *frhX*bfr&to*k<r'>!>}  - 8«C^v)l  • 

ALTROU6M  V,  - - W |T  IS  N0T  CONVENIENT  TO  EXPRESS  NETHER 

I these  relations  NOR  S-A)  IN  terms  OF  l<)  AND  14^. 

(T^ESE  ARE  EQS.  <b)  SUBSTITUTED  INTO  EOS.  aO  ( & “ ^(p  J . 


~ELE.bRAPR  EQS,  FOR  --LALA.K  t VF.LTOR  POTENTIALS  RL 


WHERE. 


<-  . C*  - . <=»*  * ^ • 

%Y  SUBSTITUTING  a-^T^Atp  INTO  THE  ALTERNATIVE  FORMS 
TCO..  ii)  CAN  ALSO  BG  WRITTEN  AS  FOLLOWS, 

,v*v  - 4j  ^ ^ 


»)  .«>  (UN|f)^(^?)  *-^^-° 

WHERE  M N.qJ.C.*  ARE  GIVEN  TO  THE  LEFT  M EOS.  F-9.). 

EQ.  it ) CAN  ALSO  BE  WRITTEN ^N  THE  FORH  : 

uy  0 * ^ fe-)  Hi  v - w 1 <£-v-° 

is  &IVSN  in  7). 


fcFQtT^ iNOEFENgENr  S^AwAN  nrrBNriAcS  V,  ^1 


ii)  (1  *N|e)V,V  - gj  ^ « E'‘)  (1  ^j^St 


I *“74 


a*»  ^.y^E  » 


tv*  as  in  ta.  7). 


'")"li  LI-.bRAPH  i a;..  FOR  THE  NOF.I'ENLEnT  sCAl-AR  pqtent'-.m  s : 

. iOC^-fc^V-^ I • 

WHERE  M,N.  cf.  c,*,  BEFORE , SiNCE  ^ • 

THE  TELE&RAPH  .TV PE  EQ.  FOR  V*  •$  ■ 

Hi)  (1  * N $'r)(  V*  IV»  - fJ^T?]  • it  S~jr. 

WWiCH  IS  NOT  AS  CONVENIENT  AS  »Q. 


llflp,^v^V^TO^TLLE^RAr»ea.  N S.ENERAL.^2l£Si^ij^~— I—i 


o>  SSS  ROE  NT, ON  y-.--s>.ro*  -fcLS»R*«‘H 


FsOLOTlON  of  EOS  9)  CO*  ABOVE  IS  WVCN  BV  , _,  . 

- (Vv)**  * r'rt 

tikew  SCALAR  FUNCTIONS  W.'*.*.  <5  * XfT,®,^,  t)  SATISFY 

0* , (I » N g? ) V*  V - a 0 *^o  h%  ■ ) 

r.  A ' . . _ A..T>  vmi  ,inr.Ai  fi\  AJiAwf  \vli  PltJD. 


THE  SO.uTiON  of  EOS.  9 CO'  IN  SShGWCAL  COORDS.  WtTU  AAIRU 
SYMMETRY  IS  NOW  Cl  VC  N BY  : V*  vt'^  • 

_ • 1 „ A M _ A LV  _ . ulv  1 


‘ £>  \{V  **J0  v',Nfe)v'3t*'  h • > 

(ij  ' OPENING  UF  T t»i  SOLUTION  LA')  AOOVE  ''V  FIND,  fjiJ 


|:?»W  * H*  sU  ^ * V 

iKAL  COA'PONENTfe  OF  T«6  VECTOR  POTENTIAL  ^ C WHICH  I S ALSO  ' 

i condition  V-  55  - o')  TO  me  a scalars  v and  x. 

Y4ERICALS  ONLY,  AHD  WHEN  IT \S  SUBSTITUTED  INTO  S«  VtpY-V.lj) 
CCMENT  - INDCPENDCNT  POTENTIAL  RELATIONS  -N  4Q4.A')  A»OvC. 

6D  TO  PROVE  this  are  : -?  V»v  - $C*V  •*)  " 

•*  ^ ^ .**■«»#»  ai*  LiiEPrtMrY  r . 4 . m • n^O  . »•*  clAEaTiCiT1* 


5UBbriru*r*NS  ^cuur^ON  (a>  int-q  ^ 

a - b f*V)]  - 11  V AS  tt,v51N  *#ovs  ’*  60&* 4)’ 


AGAIN  IN  THE  ABSENCE  OT  VISCOSITY  C1*.  W*N«*0,  «*,  SLAA j 
TlCiTYl  CO.fAl  IB  THE  SOCGNOIOAL  i^OLUTlON  OF  TVS  VECTOR. 

> 1 T J **  ^ * . *4  v .-a  _ ...  ^ ..  •vAnnt'  u<ru 


ETr  . IN  the  absence  of  viscosity  6*.  M-  N-o,  ...  sy&Tic.rO 
Binoioal  solution  or  the  VECTOR  WAVE-Ca.  |7‘-  W)|()  -o 
pIND  X SATISFY  2.  SCALAR  WAVC'SCjO  . {»«,  vy  * ^ VtO  C**o  ■••/* 


I . 1 , . - ....  - ovu^riwi^u  ttve-lp  . — 

VAVC-KQ..  (Stf  - iv  -0  IN  SPHERICAL  COORDS.  WITH 

AX.'AL  symmetry,  froV.DE 


• DEO  That  n»E  IR.AI  AR  FUNCTION 

s'at.&fTeS  A SCALAR  WAVE -EQUATION.  V'*; 

v-ors  TRAT  CO  10  *'•»')  FOR  4V  15  mCS,!**  COMPuCATSO  than 

.0,0  .HENCE,  cos.  *■>  i e>  *RE  ee&T  left  INTERNS  OF  o a u>. 


Ithe  CVL.NDRlC.AL  CASES  to  the  left  ALSO  MOLD  mere. 

QP  YHE  APARCNT  GENERALITY  OF  the  SOLENOiOAL  SOLUTION  ll)A) 

OPERATORS  * .T  'S  Only  valiO  in  spherical  coordinates,  rote 
_J  PARTS  C.«.  ARO  V»C9*C?X.)J  ) ARE  DIFFERENT  C»*, 

KenoiCular')  EVCN  when  V**-. 


lift  A PLANE- STRAIN  OR  PLANE -STRESS  SITUATION  « SPHERiCAL  COORD 
"lit  TO  THEM  is  The  AxiALLY  SvMhETRiC  CASE  0>  to  the  RIGHT. 


1 121  acvAKnU'. 

i>  REMARKS  a),  b\  «)  « r«6  CYLINDRICAL  CASES  TO  ThB  LEFT 

ALSO  HOLD  mere  . _ 

THE  MAXWELL  MODEL  >N  SPHERICAL  COORJ>S.  WILL  PC  LEFT  AS, 
u*  e»c«cise  for  the  reader. 

Hi)  this  entire  chart,  valid  for  arbitrary  t.ms-  dependence, 

CALL  BE  PAR Ti CULARllE D FOR  HARMOWC  TIME -DEPENDENCE  OP 
THE  FORH  e*U>*  IN  A COMPANION  CHART.  ft.  QAUNAUAO,  *7S).| 
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PARLE  ? POMATlONS  OF  LINEAR  DYMAM1C  V'lSCOELAS  FIG  F'Y  IN  CY=  'NDk'LAl  AND  SPHERICAL  COORDI 
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ic]  C pu>y(/»t-i<Ok„)  ■ This  ea.  CAN  6E  OPENED  UP  BY  MEANS  Of=| 

V‘tt  - [V\-  % - ^ * MeF*  - t*‘uJ  i * ' 

AND  V*.  & • 

The  gradient  operator  t the  dilatation  a are  given  above, 

ur  -9*ve*aW**?AO«£i2_]s(3.a)  +-«-!•£  «=o 

L - '<v|i,r  1 

NOTE  : V2  £ - $($•  U)  - ^nC^xU.)  . 


7)  EiElD  EOLATIONS  ; 

V1-  v(v-u)-^J)-  &iv,E.i?(vu>^i]‘|l  + 

OR  IN  COMPONENT  form,  + ia&iaA2&  U -o 

tV*Ur-^-^|^]-V^*WUr  = 0. 

HWL6  = o,  [vX]tV^Wu,.Q, 

WHERE.  Vo  . W.  wH 5a. 40-.. OMRBiw (*,*£&■) 

fjU(3«*2/i-iu>)  ’ w -irrZTZji  - <L» 

“■>■  0 - £)[  tfco.c.)  - 3*(3» so  . c.1  * & [i-oro-so- 

A - €l  'S  G'vEN  in  5s). 
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8)  Gt  RfcV.S  - NDEPENCENT  POT^NT'iAU  'TELA  T ON:- 

y)  Vl(p  + i(^-CoMir)((-f*^*-pj^<p-»^)-Kl  tr(V^). 

*u“(V  iv>^r)V^  ♦ 
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O’-*  \ ,x)1  ♦ 4^  J 

TO  WAV 

<n  mS^v,aO_-  7.'  > EOS.  FOR  r«6  scalar  l .’El-JR  PD’iNTHLC  j 

9)  MtlN'.A*  ^ l \\  1 X . i .)•<  -C  ■w.VA'ft  k . t 5k  i vlllil  M...W  . 
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WHERE,  K*  « m>*/Ca‘l1-*«>M)  . >L?  * oyc,*  (1  - IWN>  , 

NOTE  THESE  ARE  UELMhOlTJS'S  ERV  with  complex  *l‘x. 

• ) yV  * * R.‘ ! V “ o . »*  I (V*  vaJj^i-O  4. 

WnERE  , H*  » a.3.;  , -■'!».. ■ji'1  ♦-*n^iii»|AiD 

C4*  l A +H  Y-  V J 

, w* . .ifS • H-«  |Mf  sf'  v w -. 

AND,  *1  • i)!i  ♦ y.?$y  -••R./r,  vr;  - vY«.  LInLV  . 

’J>  "*r, • '••u.  # ’ o J»'  J * «,  y1.  ■_  ’ » ••  t. ‘ft  **.  ‘ l 

BBflftn 

OCVs  * . ‘O  W*  ♦ A.|  )«*>  “O  . •>•)  LVa«  *, ) JL»0  , 

WH'CRK  « IO*/iy  (!•  . M.J  « C*  Cl  — OWN) 

AND  M,N.  C,\  C,*  ARE  as  GIVEN  AB-OV6  IN  9). 

• Rf  >v  -»0  , Y»’;  < V*.  isj  ) V «V  , .I‘KV"*H,' jL“D 

WHERE  n*  AND  N,'  ARC  AS  SwEN  a&ove  N <>), 

TrE  ,‘M-T  C»«ylV  vaR  k,«y  A*>  .MlLiff.  That  r*«8  '■ATTRiAw 
IS  RVB8ER,  CVyY  *»eN  will  a.1  OCRS NC  On  ot*  A..ONE. 

“!1  ..  •■ N.  ..  * . •'  S 

t *'  *1,  ..A  v AY  ■ , . , 

"' .;  . vs 

iTlANIU  !A«UWN  *taAT  r*«E  oi.VNCIt'Ac  V <*0  > SO-vl*  sON  ->f  \V<  • aJ  )V»0  . •«*«««.»»  •»  t«H3  tOwOTiON  *?.  CJ#  i,  - *>J  w , 

ISI&IV6N8Y,  ^%a))  » (^W)*  E,  » ^vii)  - ^*".»  1 l*'  IN  CYLINDRICAL  COORDS  ONLY. 

HERE  V<*\8,  *5  ***0  Xv'  .R.iT  ARE  TWO  SCALAR  PoNCTiONN  wwtH  SATiSfV  &CALAR  '-Vv.MwCH.T a‘4,  E*&.  iO>  >•>  t i.,T. 
ta.  c*l  in  component  . form  .'ft-. 

■+V  * v |SSr  * 4*4%  • ’‘n * - ^ - ■?  liiiS  ' ■***  * |*;H  * **  £ > 

WuiCm  RELATE  tmK  TuREE  L'viRtW>CAL  C>»»C.LfeNTs  lY  rue  vL'CTOR  POTSAT**l  “i1  rC  *«.g  *-„n  scalAR  i\:nfi.;M-.  iM-  X ■ 

60eSTituTlN6  $ PUOM  f:a.(A1  INTO  11  a (?♦'  • ^.«  0 V-ELBS  TWg  OySPLACEMENT. -rDVI'CNSE **T  POTENTIAL  RELATIONS  G'yEN 
ABOVE  «H  *■>.  VOTE  THAT  V ‘S  NOT  lOt.  ;n  Tut  ABSENCE  OT  Vt&COMTV  t,« . «»n«o.  «-*.  ByARYiOTY  > Ga  Ik's  ABOVE  it 
Tnp  LO^ENOiOAl  SOlVTiOM  or  the  LECTOR  «El.Mk0lTI  Ea.  NOW  W*TR  A «?al  , RaTuKR  -HAN  COMPLEX  PROPASAViON 
CONSTANT,  MSKCE,  POR  v-8RAT,0,n's  hRmOWCiyTiMS,  S*mP  Nft  i ft  ACCOUNTED  TOR  t*v  Htg *m«,  qc  oDwPlBX  PUo'-AL AViON 
CONSTANTS  K,*  t H.J  RSLATtO  'O  T»e  El»STiC  AND  V'SCCoi  CONSTANTS  MY  ThE  POP "VL Aft  AftOvE  IN  V). 

H,  'ft  ASSOCIATED  W(Tm  hONChToOi-NAL  OR  Oil  AT  AT  O'iAL.  WAVE  S AND  M.J  Willi  IRAN  *f.R  -G  OK  t»CA«t  wAvtft. 

IN  CVi.  vt)S'C*L  COOROS.  0A&  CAN  CHECH  •rat.  V*(4«1  - AHO  "“AT.  0A  V*  v)T  a V*  v"f  . 
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13)  Nv;)  »arhG  : 

a)  ThE  KA3<t  KC.MJlT  or  *.I|J.  t.MN.r . a*0  "NAT  i>  v,hy  IT  wAS  OR-t.ANlEEO  AS  UN.*.**  *WC\€  . ;«•  That  h WL  'NTK04WC6  inDCPENDEnT  I'A'ALAK  »OVtMl 


«Ati»P(CO  « feACMONK  iH  rug  CAH-.v  vWNOt.HU',  AND  "h5N  '*<  SSST'lAUmi  NT  AND  l-ThVVS  nCvO  CO**r»NC NTS*  ARE  oiVSW  <N  UHMEi  Or  ?m£ 
PHOBCgMS.  AND  WiHCll  IS  COMMON  s.\  Ot**WK  DiSrClPGlNtii  (>«,  EGSCTN  OUvnAMu  •-  >.  ft  AHfoAHtMtCV  NOVEL  IN  NlSCOtL*&T<Ci  Ty.  AND 
W'Tu  ElECVROCSVNamU.  SkLORy  .»-ll£jl*  h\\  NOTiCE  t«£  ANAL06V  or  t«!j:  IHu  I n‘at-0M  N »mi,W  t*HlK  Wifu  TMAT  Qf  LLtCTttOMAOVUr 
ft  IT  I«  CV*t>r.NT  that  (Nall  *«*>  . A.*.J  a’*,  v‘.-**iV.  Wt  RECOVER  T^e  --imp,  ah  a1  .VLti  ir  OVNAM.V.AL  C'LAftT,ClTY  . N ThL  AUSEfiCE  < 
C)  IN  CYLINDRICAL  COORD'.-  Ill)  t . AN,  • tRt  > -.  •••  . .»  T*.  CAN  8$  <Si»C*LOr  OBTa-klO  » w.*«  tut  *’.AV6  MKAiN  At Mjl* .'•  M UHS  SJAt)  *l»OVV  0 

ti)  w twe  case  or  plane -ftT«AiN  <s&, . . a»«-  , a \ vw  o>  o a&ove)  w«ru  amal  Aymmethy  about  r«e  «.a»s,  c*«  iR-  obtanlo  »kon  su»t»  . 

«)  A CASE  or  SO*  .NTERESV  vyu^h  REt'Ov  t Twt  Iwv  ( sCv'oS  CONSTANTS  TOONE  .OCCURS  V4L  "«E  ttCH.-V.SCOt.iTY  < EU.O  OR  Hi  ALLY  SO.  "-  i AM 
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j^4»  - r a*  - v_ai  icj  *■? 
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^Ii2  v?) _ ^ ^ ‘ ** 

TO  HAVE  THCSf  IN  TYRkS  OF  Cf  4 !•>,  WE  REPS-ACE  ■'J1  W V> 


I A,-  -•  -.  • '.«i+  - *.*•[&  • *1  i>(Val)]]  + ^ 

1 M,i-A.N^  ♦***[?&♦*$&  - k*MtH)]}  +m> 

j tjj”  M«  {-^«K.f<p}  ♦ N,  , t ®«r  “° 
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i H sO'-rAST  1 '..=  REPLACE  :<.iT-  8V  ’,)* 
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r.f  - ■4>i/»-,V*»‘^M') 


- K.1 ',  'P  « O 
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i a ~ -Me  y ' 

«£*E - -£.V*X(r,e)  - WGNCC  U)  ACOVE 

»5  EQ'Jl  valfnt  to  cither  ALTERNATIVE  form  fc*VEN  IN  io). 


4 2 » ^xi  ^ 4 -» ( 3«  ♦ 4-tO  -3/  — _ V.  4. 

— m4  T ~ ~ J .*2 1 /»*//>*  ^ a _ 


;0  (v*+k?)V  -O  l v-Ip-o 

ilzsiu 

J*  J . _n-a/-3  3 a _ 


« ; (V  * ♦ H.J  ? <p  - O t'i  ) (v2  * j "+4  “ O 

WHERE  , >Sj‘,  M,  N,  C/,  c.*,  ARE  AS  GIVEN  AbOVE  M O). 
,-ts  is  EQUIVALENT  TO  ! _ ^ [(^+  ^ mQ  , 


o (v2  ♦ k;4;<p  = o ;o  (v1+kJ)  v,*o 

WHERE  rtf  l H-i  ARE  AS  GIVEN  ABOVE  IN  9).  SINCE  V,-  K^JS 
Ea.  *0  is  equivalent  to  - v*  [(v*t  -o. 

the  CONDITION  C,*«Ci  MEANS  The  MATERIAL  16  RUBBER. 
Only  for  c»«cj  will  K,1  DEPEND  ON  C 1 ALONE. 


■n  ^ulcnoical.  ;jmriO.N  c-'  \.,  _r-jR  l^nu-o'jtz-s  equation  >n  --anc  ivun  coo^uinvls  : 

-HE  iD-LlTiON  -T-'lO'wi  -WAT  SiVEN  TO  THE  -EFT  FOR  oEKLKAL  n'UNCRiLAL  COORD.nATES  BUT  NOW  WITHOUT  THE  f DEPE N DENC 1 
hence,  one  can  also  show  that  -he  solenoidal.  solution  of  CvS-«4^*  =0  ! which  is  also  the  solution  of  “*C”W)- 
-K^V  .’  IS  GIVEN  BY,  s . <L«  e.Ki  X O) 


-H^V  . 'S  GIVEN  BY,  Vp  - ^r.«  (Vr,*(2.X^-  e,K i X O) 

where  xCr, e)  is  a scalar  function  satisfying  the  scalar  he lm-oltz.  equation  (t'+kIOX-o.  equation  Ca~)  in 
component  FORM  is  . V,-0  . V.-0  , l«,  WN'Ch^  relate  the  cylindrical uc«»NBNrs  OF  lJ 


to -he  scalar  function  X-  Substituting  solution  (a)  into  Tt  = 9<p*VxJ  yields  a =^<p-  h^(5?x)ke»  which  are  the 
C SPLACFMENT-  independent  potential  relations  found  in  eqs.  a")  above. 

UE  REPEAT  THAT  IN  the  ABSENCE  OF  VISCOSITY  0«-,  M-N-O  , vi»,  ELASTICITY),  SOLUTION  (AT  IS  THEN  THE  SOLENO'DAL  SOLUTION 
OF  THE  VECTOR.  HELMHOLTZ.  EQUATION  WITH  A REAL  RATHER  THAN  COMPLEX  PROPAGATION  CONSTANT,  PROVIDED  THAT  X SATISFlC 
A SCALAR  HELMHOLTZ  Eft.  WITH  THE  SAME  REAL  PROPAGATION  CONSTANT"  4L*  - W/c4  . HENCE  , DAMPING  DUE  TO  VISCOSITY  1$ 
ACCOUNTED  FOR  IN  THESE  MODELS  BY  MEANS  OF  COMPLEX  WAVENUMBERS  X*  and  K.J  RELATED  TO  THE  MATERIAL  PROPERTIES  AS 
SHOWN  ABOVE  IN  9).  K,  U:  ARE  ASSOCIATED  WITH  LONGlTUOINAL  AND  TRANSVERSE  WAVES  RESPECTIVELY.  NOTE  THAT  in  THE 
MAXWELL  MODEL  THE  PREVIOUS  STATEMENT  ‘S  TRUE  ONLY  WHEN  C{«c£  . OTHERWISE  K, "CONTAINS  BOTH  TYPES  OF  WAVES. 


AR  POTENTIAL'S  WHICH  SATISFY  The  HELMHOLTZ'S  ECJ.G.  WITH  COMPLEX  *.'»  IN  IO) , THEN  THE  FIELD  EQS.  7)  ARC  AUTOMATICALLY 
IMS  OF  the  POTENTIALS  SY  E 05.  4)  A 6)  RESPECTIVELY  , THIS  TECHNIQUE , WHICH  IS  BASIC  T0  SET-UP  AND  SOLV&  BOUNDARY  VALUE 
CITY.  AND  THAT  IS  WHY  ~HlS  CHART  WAS  DEVELOPPED  FOR  ALL  THESE  CASES  AND  VISCOELASTIC  MODELS.  HE.  READER  FAMILIAR 

romagnetic  theory  in  a coulomb  rather,  than  a lorentz  (<t,  v-  4)  ♦ £■  f£  gauge.  > 

ABSENCE  OF  Viscosity  (*«.  A„.-  JU„»o  =or  the  kelvin-voigt  model  and/or  a»p-o  for  the_  maxwell  model .) 

^ AIJOVE  UY  the  STANDARD  TRICK  of  REPLACING  X,  IN  THE  PLANE -STRAIN  CAST  I.  BY  A,  DEFINED  TO  BE  I A*  - !AA«/WC-'*+  2/*0  • 

M DUE)  BY  SETTING  W..O  AND  **  CANY  VARIABLE  J • O , THEN  THEKF.  IS  ONLY  r- DEPENDENCE  . Ci«,  Ur*  ILtr.)  8 Kf-<pCr)  l .Ur*  W/W  ■ 

JO  -HIS  AMOUNTS  TO  HAVING  3A^r2/4B-0  (HKLVIN- VOIGT)  OR  Oor*  2p  *0  (MAXWELL).  SO  THE  TO^MULAft  THEN  SIMPLIFY  A LITTLE. 
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WE  WKL  SEE  BELOW  ThAT  IT  ,S  not  wO 
terms  O f if  i v,  - • •>.£  .‘MCE  /p.  LAI 
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""6c:e  »Kc  nEUMmOcTZ  EOlS.  with  COMPLEX  PROPAGATION  CONSTANTS.  The  vector  EQ..  IS 

plrtmer  Subjected  to  the  solenoidal  Cgauge)  condition  shown. 


„,.  - A ' "VNbc  ■ N Jt  *-  Ci-iPtN  r v-v.ALA^  Tfa  .'4T  i/a  LS  • 

) (V*’  - K,a)(p  = O , CVJ"K|)!p=0  , CV2♦'<J)X^.r.e.*■)  -o  . 

WHERE  The  COMPLEX  -v1  & "-1  *PE  AS  Given  IN  9)  AND  VJ  AS  GIVEN  IN  7) . 
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CT  (7-  . W-  ) V - 9 j "«C  ^WUTiON  Of  v«tf»  ip)  - KJ  ip  ANDiTIS, 

« Cjrl,  r*  V H Curl cwrL  * t r - vx  - K?  r X * C (VZ.  • ?)  -^V>*T  + ^ (A> 

PROVIDED  Ipi.'".  0.<»  i Xt ",  0,  SATIDPY  SCALAR  HELMHOLTZ  EQ.S.  (V1*  K|  ) !p  C«r  X)  =■  O . 

SOTE  "HAT  Ip  IS  NOT  i V I • OPENING  SOLUTION  l>}  INTO  ITS  COMPONENTS  WE  FIND, 

viCTZ)*  r«.*X  ■-lii.vala  it  . ,_i_  ^ “ stW  ^ 4 V Si  + Ms  • fBs 

sr*  i-35i  r W P ,iu  TJT  ’ M . 3^  _ ! i*  3aX 

which  RELATE  Vr.4>»,  TO  Ip  AND  ” * 3e  i"i.ne  Tir  bine  or*>  ’ 

J5  NG  THE  AUXILIARY  RELATIONS  V-  O^V)  ] * V<(J  » "j  ? V " 9 CC  4.'  • r)  AND  V «.  (P  X.)  -(VRC)*  r WE  CAN 

now  Substitute  solution  (a)  into  U » Pip  t 9« -p  to  obtain  u . i?  [<*>+  l^ov)]  + Kj4  [ri^i  ♦ 
as  we  mad  1 n eo.s.4.  in  the  absence  of  viscosity  0«,  m-n-o  , elasticity)  ea.CA)  above 
is  "he  SOLENOIDAL  SOLUTION  of  the  vector  HELMHOLTZ  E.Q.  NOW  WITH  AREAL  RATHER  THAN 
COMPLEX  PROPAGATION  CONSTANT  Af  * ■O’/cf  PROVIDED  Ip  A X SATISFY  SCALAR  HELMHOLTZ  EfiLS 
WITH  THE  SAME  REAL  PROPAGATION  CONSTANT  *1  . 
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these  are  helmhol-z  eas.  with  com 
since.  <p«.  u\^.d)it  and.  v^(ip4^) 
turns 'our  that  the  eo  tor  Is  ••  ^ 
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WHERE  Kr  L k|  are  The  COMPLEX  QUAN 
IF  WE  SET  - - SV  m the  uAfT  EG.  Gl 
This  IS  WHY  IT  is  MUh»  CONVE  NIENT  to  i 
terms  of  p rather  than  4). . 
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the  solenwidal  w-.Lvr'ioii  jp  (vr;*"j 

if)  * ( vrjl  4) ) a i*  ■•=  R?  r (P  -J(r,  6)) 

OR  IN  COMPONENT  FORM.  VV“°  • V# 
( NOTE  That  :p  16  NOT  ' p|  ,)  PROVIDED  "ha; 
SUBSTITUTING  L-OLUTi  ON  CA  ) INTO  Irn 
ft  * ’■  Jr  Vr  4 .:t  Ip  AV-  WE  HAD  A 

SOLUTION  r A)  MOLDS  '\  ThE  ABSENCE  of \ 

elasticity''.  ,t  is  then ~i-e  solution  o> 
E&.  with  a heal  Propagation  conctai 
SATISFIES  THE  SCALAR  HELMHOLTZ  EG 

gation  constant,  note  *»at  the  la 

HELMHOLTZ  EG.  WITH  KJ  - K,*  - 1/r-».Si 
TO  SOLVE  THAN  IC'Ki'}  BUT  IT  IS  L-TILl 

IS  possible  FOR.  «.*■  AS  general  AS,  V 
* rlh?3*W)  WHERE  f.4,  K,  ARE  CONSTA 
FUNCTIONS  Or  71 1 K INDICATED 


•■il  uLMAHj  ■ 

.>  RL marks  a.),  b),  «)  OF  THE  CYLINDRICAL  cases  TO  THE  LEFT,  ALSO  HOLD  HERE. 

u)  NOTE  AGAIN  That  IN  SPiTE  of -he  APPARENT  GENERALITY  OF  THE  SOLENOIDAL  SOLUTION  1l)A>  ABOVE 
Wi?M  ITS  VECTOR  OPERATORS,  IT  Ifl  ONLY  VALID  IN  SPHERICAL  COORDINATED.  NOTE  FL'RTHCR.  THAT 
the  TWO  PARTS  C'««  Cnrt  I*V  AWD  Curl  Curl  i1*)  ARE  DIFFERENT  EVEN  WHEN  V-X.  ACTUALLY, 
They  ARE  PERPENDICULAR.  WHEN  V*X.  , 

;.<) There  s noway  to  have  a plane-strain  on  plane-stress  situation  in  spherical  coords. 

THE  CLOSEST  wt  CAN  GET  TO  THEM  IS  THE  AXIALLY-SYMMETKIC  CASE  0)  TO  THU  RIGHT. 
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C)  REMARKS  a),  t)  t «)  OF  THE  CYLINDRI 

also  hold  here,  remark  a>  is  t»e  * 
vi)  The  maxwell  MODEL  N spherical  CC 

A9  AN  EXERCISE  TO  THE  READER. 

A'O  this  entire  Chart,  which  >s  valid 
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. TmAT  .-0  -«E  “AWtLL  N ^IMLR-.-AL  COORDINATES  Wii_L  PC  LEFT 

JALlY,  AS  AN  LXEHL I t E tO  Tug  READER, 

Wi;  Tm!Ji  entire  Chart,  ■„■.<, S VA_,D  rDR  mAN*VJ\*C  time -DCREN- 
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r ?•-  4»-  Jf*.  “men  -MEkU  ARCTVO  ALT^RNATr.C*, 
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" ' ^iir(.  lir  y*  j l-TSneJ  * ^ ( K-vi- >J  1 11”  V*  > ' 404  n®^ 

FOR  S>0  OR  5*0  RESPECTIVE WY.  >NTEC,ER,“)  . ^ 


. • ,is.  c<  : -■,.-..OOC»M  T-t  -CcMNOOE  OF  INTAODOClNi  •NOEReNDeNT 
^LA.AR  Rjri'.r.A.S  WM.CM  j»T|spr  TELECOAPh  Eas,  Cor  -ELMMOLTl'S 
CDS  ViTM  CC^flEY  -S30FA&ATIDN  CONSTANTS.)  SO  thAT  ALL  TME  STRES- 
SES l D-SPLAtSMEN'e  ARE  CCTERMiNEQ  ATIOM  TmOSE  POTENTIALS., 
■'■  APPARENT lY  novel  IN  VlSCOKLAf-TiCiTY,  -MERE  'S  much  OTHER. 

. -c mature  >1  VSCOSlASTil  models.  «^or  example, 

-)  W TOLD  NOVA.*  1.'  DYNAMOS  Of  ElASTC  SYSTEMS'.'  TOHN 
U-EV  i SONS.  INC,  I RfcO  CCHAR  2) 

_>  A CEMAL  CRINSEN,  - MECHANICS  OF  CONTINUA"  ODMN  VIILEY  S 
SONS.  INC  . ’9ST  , CCHAPTEK.  9). 

VCT  Mix,**  acpl’EO  -ATERIAL  is  &IVEN  IN  THESE  REFERENCES 
?ETO*JD  Tm£  field  EQUATIONS  t CONSTI TUTI VE  RELATIONS, 
JSmA-.lv  m GENERAL  TENSOR.  FORM.  OTHER  MORE  toENCRAL 
V‘Co'ELAt-C  "ODELS  CSUCm  AS  THE  'STANDARD'  VI5CO  - 
Elastic  y o del  .which  is  a parallel  combination  of  a 
ka/wELL  and  A KELVIN -VOKST  MODELS  AT  EACH  MATERIAL  > 
POSIT  OF  THE  BODY,  AND  WHICH  IS  VERY  COMPLICATED  TO  HANDLE) 
ARE  O'SCUSSED  IN  6ENERAL  'N  ThE  ABOVE  REFERENCES. 

SOME  OF  TmE  BASIC  EOS.  OF  STATIC  ELASTICITY  C!».  TlM®* 
DiPENOENCE  t NO  VISCOSITY ) N £|K  ORTHOGONAL  COORDI- 
NATE SYSTEM-  mane  SEEN  DEWIVLO  Py  STANLEY  HALPER60N 
M S.TmE.miS,  CATHOLIC  UNIVERSITY,  WASH.  D.L.  , APRIL  I9ST  .)  . I 
THIS  ISONC  OF  THE  NICEST  COMPILATIONS  OF  THIS  .jORT  WG  j 
HAVE  SEEN. 
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in  a more  or  less  complicated  fashion  depending  on  the  model  being  used.  For 
example,  the  relation  is  more  complicated  for  the  Maxwell  than  for  the  Kelvin 
model,  in  any  coordinate  system.  These  tables  provide  us  with  the  methodology 
and  information  needed  to  solve  boundary-value-problems  of  viscoelasto-dynamics 
for  spherical  or  cylindrical  geometries,  in  the  Kelvin-Voigt  or  Maxwell  models. 
The  approach  and  material  covered  in  these  tables  seems  to  be  novel  as  applied  to 
viscoelasticity.  The  reader  familiar  with  electro-magnetic  theory  will  quickly 
notice  the  analogy  of  the  presentation  in^this  table  to  that  of  electrodynamics 
theory  in  the  less  familiar  Coulomb  (v  • = 0),  rather  than  Lorentz  (ie,  V • \p 

+ c It  = 0)  gauge* 

To  verify  that  all  the  entries  in  the  tables  check  with  each  other,  we  will 
illustrate  the  use  of  the  formulas  by  analyzing  one  example.  We  arbitrarily 
select  column  D,  which  deals  with  the  Kelvin-Voigt  model  in  cylindrical  coordi- 
nates for  the  subcase  of  ^ .ane-strain. 

1.  Definition  o..  Plane-Strain 


In  piane-str- in,  the  strain  tensor  is  two-dimensional  and  all  strain  compo- 
nents with  a subii.dex  z (axial  coordinate)  vanish.  The  stress-tensor  is  not,  and 
there  is  a nonvanishing  normal  stress  T . Another  equivalent  way  to  define  this 
situation  is  by  stating  that  the  axial  SJsplacement  u vanishes,  and  that  the 
other  two  displacement  components  uy  and  u^  are  z-independent.  Plane-strain  is  a 
particularly  useful  approximation  for  bodies  which  are  very  long  in  the  z-direc- 
tion. 


2,  Constitutive  Relations 


These  are  the  relations  between  the  stress  and  strain  field  components. 

There  are  three  normal  and  one  shear-stre:..'  components  related  to  two  normal  and 
one  shear-strain  components.  Stress  and  strain  are  independent  concepts.  Only 
when  these  two  notions  are  linked  through  constitutive  relations,  is  that  a visco- 
elastic theory  is  formed.  The  first  of  these  four  relations  is, 


T 


rr 


= ( 


X 

e 


iu)Xy)  A + ?(ye  - iwyv)  err 


where  A is  the  trace  cf  the  strain  tensor  (ie,  £ rr  +epQ^  ^ is  ca^e(i  the 

dilatation.  The  other  three  relations  are  shown  in  Table  2. 


3.  Strain-Displacement  Relations 


Since  the  elastic  or  viscous  constants  do  not  appear  in  these  relations, 
they  are  the  same  as  in  elasticity,  namely, 

3u 

Y 

e 


rr 


3r 


e06  55  i + !x  , 2e  « 1 9ur  » 3u6  - u6 
0 r 36  r r8  r 36  3r  r 


all  other  strain-components  vanishing  by  definition  of  plane-strain. 
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h . Displacement  - Independent  Potential  Relation 

These  relations  are  intimately  linked  to  the  solenoidal  solution  of  the 
vector  Helmholtz  equation  in  plane  polar  coordinates  which  is  discussed  below  in 
item  (ll).  The  relations  can  be  written  as  a single  vector  equation  as  follows, 

u * ^ ) x t 

z z 

where  <J>  and  ip,  are  the  two  independent  scalar  potentials^needed  in  this  plane- 
strain  formulation.  It  is  shown  in  item  (11)  that  t|>  - x>  vhere  x is  a 

solution  of  the  scalar  Helmholtz's  equation  ( ) X * 0. 

5.  Stress  - Displacement  Relations 

If  the  strain-displacement  relations  (3)  are  substituted  into  the  constitu- 
tive relations  (2),  the  result  is  the  tress-displacement  relations.  The  first 
such  relation  is, 


«■  ^ 

3u 

T = [a  - iiuX  | 

A + 2 fy  - 

itoy  ] 

r 

rr  1 a v] 

Le 

VJ 

3r 

where  the  dilatation  A is  now  expressed  in 

terms 

of  the  displacement  component 

as  follows, 

3u 

, 3up. 

u 

A - 

+ 1 0 + 

r 

• 

* “ 3r 

r 38 

r 

6.  Strain  - Independent  Potential  Relations 

If  the  displacement-independent  potential  relations  (1*)  are  substituted  into 
the  strain-displacement  relations  (3)»  the  result  is  the  strain-potential  relation. 
The  first  of  these  three  is, 


and  analogously  for  the  other  two.  Here  $ and  ^ are  the  two  independent 

scalar  potentials.  Not^  that  \p9  ^is  the  axial  component  of  the  vector  potential 
Further,  ij>  = X 55  X vhere  X is  another  independent  scalar  poten- 

tial vhich  can  also  be“used  here.  (See  item  (11)  below) 

7.  Field  Squat  ions 

These  are  the  Navier  equations  given  initially  it  (i).  For  harmonic  time- 
dependence  they  can  ho  written  as  follows, 
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T u + 


1 + 


A -iaiA 
e 


-iaiA  I 

-r-i  ? ft  • Z 

-io )y 

V J 


) + k2u  3 0 


2 2 X ~fc 

where  <„  = pco  [ y . - iuy  ] , and  A = dilatation  = V • u.  For  plane- 

strain  in  cylindricals , the  quantities  V u and  V A take  on  the  following 
simplified  forms, 

= u - _r  - d el 

v r 2 2 ae  J 


v2  u 


£ + 
r 


{v2ue  - !i  * 2_  !!tl 
10  o 2 ae  ; 


, a 3A  ft  ^ 1 3A  g. 

grad  A ” r-  e + — • r^-  6 

3r  r r 30 


and  also  we  have. 


u = ur  Er  + u8  H 


and, 


/a2  + i a_  + jl  a \ u (r,e). 

r \3r2  r 3r  r2  382J 


8.  Stress  - Independent  Potential  Relations 


If  the  strain-potential  relations  (6)  are  substituted  into  the  constitutive 
relations  (2),  we  obtain  the  stress-potential  relations.  One  form  of  the  first  of 
these  relations  is, 

V = [xe  ' tuXv]  * * 2 [ue  - luuv]  {^£  * If  [r  5^]} 

where,  A 3 * u 3 $ • x I];]  = ® -k:^2  <£> 

2 

The  second  form  can  be  found  by  setting,  >i>  ° k.  y . The  formulas  for  the 

Z w 

other  three  stresses  present  in  this  case  can  be  found  in  Table  2. 

9.  Helnholta's  Equations  For  The  Scalar  And  Vector  Potentials 

Substituting  u 3 grad  $ + curl  into  the  Navier  equations  vritten  in 
the  alte.*  i‘ < *ve  form, 

{X*  ♦ 2y*)  $ • u)  - u*  s (§  x u)  + p*  u 3 0 

£. 

2 2 

where  X*  = X - iuX  , y*  3 u - iwy  , and  w*  <_  3 pu  we  eventually  find, 

6 v C V c 
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(A*  + 2y*)  + k:^2^}  + y*  x {v2^-  (§  • "^)  + 11)}  = 0 where 

2 2 2 
k,  = pw  [A  + 2p  - iw  ( A +2y  )]  ~ and  K.  is  as  we  defined  it  above  in 
1 L e e v v J 2 

item  (7).  It  is  obvious  that  the  above  expression  is  satisfied  if, 

V2  4 + = 0 and  V2^  + K^2  4)  = 0 provided  that 

• Ip  = 0 (ie,  "Coulomb,t  gauge). 

10.  Helmholtz's  Equations  For  The  Independent  Scalar  Potentials 

p p 

Since  in  this  case,  \p  - £ tj)  then,  V \p  = £ V \p  , and  the  result  is, 

Z Z 2 2 

(v2  + Kx2)  4)  = 0 , (V2  + <22)  >Pz  = 0 

11.  Solenoidal  Solution  of  The  Vector  Helmholtz  Equation 

2 2 -♦* 

It  is  not  hard  to  show  that  the  solenoidal  solution  of  (V*~  + ) ij)  = 0 

is  given  by,  Ip  = ^ x x (6  y)]  = Ko^  X where  x (r,0)  is  a scalar 

2 z c g g 

function  satisfying  the  scalar  Helmholtz's  equation  (v  + k.  ) x = 0.  In  compo- 

c p 

ent  form,  the  above  relation  is  ^ =0,  \p„  = 0,  and  i|)  = kl  x • These  are  the 

relations  between  the  three  cylindrical  components  of  t£e  vector  potential  ip  and 
the  scalar  function  x (r,8).  Note  that,  indeed,  v • ip  = 0.  Note  that  in  the 
absence  of  viscosity  (ie,  A^  =»  0,  yy  * 0)  the  complex  propagation  constants  ic ^ 
and  become  real.  Thus  viscous  damping  is  accounted  for  in  these  models  by 
complex  propagation  constants,  related  to  the  material  "constants"  as  shown  above 
in  items  (7)  and  (9). 

12.  Remarks 


a)  In  cylindrical  coordinates  the  plane-stress  results  are  derivable  from 

the  plane-strain  results  by  replacing  the  elastic  constant  A^  in  the  plane-strain 

results,  by  the  fictitious  "constant"  A defined  to  be,  A = 2A  y j(A  + 2y  ) . 

e e e e e e 

In  spherical  coordinates  there  is  no  way  to  define  the  plane-stress  subcase.  In 
plane-stress,  the  stress  tensor  is  two-dimensional,  which  means  that  all  stress- 
components  with  a subindex  s,  vanish.  The  strain-tensor  is  not  two-dimensional 
and  there  is  an  e nonvanishing  strain.  This  case  is  ideally  suited  for  bodies 
which  are  very  thfH  in  one  direction,  the  z-direction. 

b)  Plane-strain  (or  plane-stress)  with  axial-symmetry  about  the  z-directior. 
can  be  obtained  from  the  results  in  the  eleven  items  above  by  merely  setting  u^  = 
0 and  3__  [any  variable]  » 0.  There  is  only  radial  dependence  in  this  situation. 

38 
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c)  All  the  above  applies  to  the  Kelvin-Voigt  model  in  cylindrical  (actually, 
plane  polar)  coordinates,  for  the  plane-strain  subcase.  This  is  Column  D of  Table 
2.  There  are  twelve  cases  covered  in  Tables  1 and  2 and  this  one  was  intended  as 
an  example  to  show  how  the  various  entries  are  derived  from  the  others , how  they 
check  with  each  other,  and  how  one  could  proceed  to  drive  any  other  case.  We  have 
used  the  information  in  these  tables  to  rigorously  set  up  and  solve  acoustic  scat- 
tering problems  where  the  scattering  bodies  are  elastic  cylinders  and  spheres  with 
their  outer  surface  coated  with  layers  of  viscoelastically  absorbing  materials.4 
One  reason  to  have  written  this  report  is  to  gather  these  fundamental  viscoelastic- 
ity relations,  which  are  not  available  elsewhere  to  this  degree  of  detail  and 
approach,  in  one  document  that  we  could  refer  to  in  future  work  as  "he  place  where 
the  theoretical  background  is  derived  and  presented,  in  the  form  in  which  we  will 
use  it. 

d)  Some  authors  have  stated  that  in  order  to  solve  problems  involving  absorp- 
tive bodies  one  must  deal  with  Helmholtz's  equation  with  a complex  propagation 
constant  k.  This  is  indeed  correct  but  there  is  much  more  to  it  than  just  that. 
Just  with  a complex  < we  would  not  know  how  the  real  and  imaginary  parts  of  K are 
related  to  the  elastic  and  viscous  material  constants  of  the  solid  in  the  various 
viscoelastic  models  that  one  could  use.  In  fact,  we  would  not  know  this  relation- 
ship in  any  model.  Furthermore,  we  would  not  know  how  to  relate  the  stresses  and 
displacements  in  the  body  to  the  solution  of  that  Helmholtz  equation  with  a com- 
plex k.  Hence,  although  the  idea  is  correct,  in  practice  one  really  needs  to  de- 
rive all  the  detailed  information  contained  in  these  tables,  and  that  is  why  we 
developed  them.  Careful  examination  of  these  tables  shows  that  the  problem  is 
really  harder  than  anticipated,  since  we  must  solve  not  one  but  several  (ie,  three) 
Helmholtz’s  equations  with  various  (ie,  two)  complex  propagation  constants  which 
are  different,  and  then  go  through  various  other  sets  of  equations  (ie,  1 and  8) 

to  obtain  the  displacements  and  stresses  from  the  solutions  of  those  Helmholtz's 
equations.  We  finally  point  out  that  this  procedure  yields  different  results  in 
each  one  of  the  various  models  and  cases  presented  there. 

e)  Column  A,  Table  1 shows  some  general  formulas  valid- for  all  coordinate 
systems.  Note,  however,  that  since  the  solution  of  the  vector  telegraph  equation 
varies  with  the  coordinate  system,  not  many  general  entries  can  be  filled  in  that 
column. 


* G.  C.  Gaur.aurd,  "Sound  Scattering  from  an  Elastic  Cylinder  Covered  With  a 
Viscoelastic  Coating",  JASA  58,  S1C1,  1975.  Also,  Free,  of  III,  U.S.-Fed. 
Republic  of  Germany,  Hydroacoustics  Symposium,  Munich,  Germany,  Vol  1,  Fart  II, 
PP  UU,  May  1975  (U) 
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or  also,  <Mx,t)  = A KX+i0Jt 


where,  < = - 


Cp  (1  + iuP) 


Here  C is  real  and  the  - signs  in  front  of  ioc  describe  waves  travelling  to  the 
right  or  left  of  some  origin,  (if  one  were  interested  in  time-dependence  exp(-iwt) 
we  would  replace  i by  -i  in  all  these  results.) 


Now  let  us  consider  the  standard  one-dimensional  wave  equation, 


2 2 

a 4>  * 1 31 

2 2 o 

ax  c at  • 


The  solution  of  eq.  (4)  for  harmonic  time-dependence  exp(iojt)  can  analogously 
be  written  as. 


4>(x,t)  •*♦<*]  where  <2  . S- 


or  also. 


+ kx  + iwt  . 2 ui 

(p(x,t)  - A e-  where  < = -or 


It  is  obviouSgthatjSOlutions  (2)  and  (5)  and  also  (3)  and  (6)  can  be  made 
equal  by  setting  (l  + iwP).  It  is  clear  that  we  can  solve  the  telegraph 

equation  with  a real  propagation  speed  C by  ignoring  the  damping  term,  ie,  by 
solving  the  wave-equation  with  a suitably  chosen  complex  propagation  speed.  Com- 
plex speeds  are  artificially  produced  by  fictitious  complex  elastic  constants  that 
one  can  introduce  for  this  purpose. 


Let  us  now  look  at  the  case  of  shear  waves  in  a solid.  Define  a complex 
shear  modulus  u*  ■ U*  ♦ iu”  * u*  (1  ♦ 14),  where  6 » u'Vu*. 


The  equations  of  linear  viscoelasticity  are  the  same  of  those  of  elasticity 
(ie,  no  viscosity)  if  X^  and  in  the  elasticity  equations  are  replaced  by  the 


quantities,  X*  » X + iwX 

e v 


u*  ° u ♦ im 

e v 


NOTE:  This  is  true  only  in  the  Kelvin-Voigt  viscoelastic  model  with  assumed  time- 

dependence  of  the  fora  exp(iut). 


Clearly  U = U‘  and  u a u"/a>.  Thus,  for  shear  waves, 
e v * 


p B 3 a 7T  3 frr.  = “ and  C - C /TV 
U U)U  U)  s 


V ± lu" 

p 


VF 
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Hence,  for  harmonic  shear  waves  in  the  Kelvin  model,  by  solving  the  wave- 
equation,  with  a complex  shear  modulus,  we  have  the  solution  of  the  telegraph 
equation.  This  can  also  be  done  in  the  Kelvin  model  for  longitudinal  (or  dilata- 
tions!) waves. 


It  turns  out  that  the  Kelvin  model  with  harmonic  time-dependence  assumed,  is 
the  only  case  where  the  viscosity  coefficients  X , y^  of  the  viscoelastic  solid  are 
linearly  proportional  to  the  imaginary  parts  X",  y"  of  the  "complex  elastic  con- 
stants" X*,  y*,  the  proportioned ity  factors  being  1/w.  It  is  also  the  only  case 
where  the  elastic  constants  X , y of  the  solid  are  just  the  real  parts  X',  y*  of 
the  "complex  elastic  constants." 


To  show  this,  let  us  now  consider  the  Maxwell  model,  where  the  springs  and 
dashpots  at  each  material  point  are  now  connected  in  series.  It  can  be  quickly 
shown  that  the  field  equations  and  constitutive  relations  of  the  Maxwell  solid 
are  the  same  as  those  of  elasticity,  provided  that  the  elastic  constants  X , y^ 
of  the  elasticity  equations  are  replaced  by  the  operators, 


X , It  jXe(2B  * lt°~  2ctye  ) 

26+  It  * 301  + 26  + It  ^ 


3 

It 


26+lt 


(9) 


where  a,  6 are  the  viscosity  coefficients  of  the  Maxwell  model  (which  are  analo- 
gous to  X , y of  the  Kelvin  model).  For  time-dependence  of  thi  fora  exp(-iwt) 
we  can  call  tKese  operators  by  the  name  "complex  elastic  constants"  X*,  u*  ie. 


X* 


IB  (V29  - lM>  - ^ 
iw-28  I 3a  ♦ 28  - iw 


U* 


yftiu) 
its  - 28 


(10) 


It  is  obvious  that  if  there  is  no  viscosity  (ie,  a = 0,8  = 0)  these  quantities 
X*,  y*  would  reduce  to  X , u respectively.  It  is  also  clear  that  a and  8 are 
now  not  proportional  to  tRe  imaginary  parts  of  A*  and  y*  respectively.  Further- 
more, X and  u are  not  the  real  parts  of  X*  and  y*  anymore.  Thus  the  "trick"  of 
the  complex  elastic  constants  does  not  work  here  at  all*  The  ^rick*  also  fails 
for  the  standard  viscoelastic  model,  or  any  more  complicated  model  which  contains 
at  least  one  Maxwell  element.  In  fact  it  fails  even  when  there  is  no  Maxwell 
element  in  the  model  provided  there  is  more  than  one  Kelvin  element.  Thus,  only 
for  one  single  Kelvin  element  will  it  "vork." 


For  the  reasons  given  above  we  believe  that  the  most  clear,  systematic  and 
natural  way  to  handle  viscoelastic  problems  of  this  sort  is  to  measure  the  (real) 
elastic  constants  of  the  material  Independently  from  the  (real)  viscous  coeffi- 
cients and  then  use  those  numerical  values  in  the  field-equations  when  analyti- 
cally solving  them.  Of  course,  ve  must  keep  the  damping  terms  in  the  field  equa- 
tions. The  final  result  will  be  complex,  but  this  fact  is  due  to  the  presence  of 
the  damping  terms  rather  than  because  of  any  "complex  moduli"  we  want  to  ficti- 
tiously introduce  because  of  our  insistence  on  ignoring  the  damping  terms  instead. 

We  have  seer,  that  the  complex  moduli  "trick"  becomes  exceedingly  difficult,  in  fact 
impossible,  in  any  case  other  than  the  Kelvin- Voigt  model  vith  harmonic  time-depend- 
ence. 
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1.  The  Correspondence  Theorem 

A very  useful  theorem  that  permits  us  to  relate  solutions  of  elastic ^ problems 
to  those  of  viscoelasticity  (in  any  model)  is  the  Correspondence  Theorem.'  This 
theorem  has  a conceptually  very  simple  statement,  hut  it  is  very  difficult  to 
apply  it  to  actual  cases.  The  theorem  basically  says  that  if  we  want  the  solution 
to  a viscoelasticity  problem  (static  or  dynamic,  and  in  any  model)  we  should  first 
solve  the  "corresponding"  elastic  problem  without  viscosity.  Then  we  Laplace- 
transform  the  solution.  Then  we  replace  the  elastic  constants  appearing  in  that 
Laplace-transformed  solution,  by  certain  "memory  functions"  (which  vary  with  the 
viscoelastic  model  being  used).  The  result  of  this  replacement  is  the  Laplace 
transform  of  the  solution  of  the  "corresponding"  viscoelastic  problem.  Inverting 
it,  we  have  the  solution  of  the  viscoelastic  problem  we  wanted  to  solve.  We  here 
recall  the  Laplace  Transform  pair, 

F(s)  = £°°  f(t)  e”stdt  , f(t)  = ^ i“  F(s)eStds  . 

J cr-i00 

The  prescription  given  by  the  Correspondence  Theorem  is  very  straight  forward. 
It  turns  out  in  practice  that  when  the  elastic  constants  are  replaced  by  those 
"memory  functions",  the  resulting  Laplace-transformed  solution  that  must  now  be 
inverted  is  quite  formidable  in  most  cases  of  interest  (ie,  dynamic  cases).  We 
should  point  out  that  the  memory  functions  depend  on  the  Laplace  transform  variable 
s in  a more  or  less  complicated  manner  depending  on  whatever  viscoelastic  model  one 
uses.  See  equations  (7)  or  (10)  with  iw  replaced  by  s,  for  the  Kelvin  or  the 
Maxwell  model  respectively. 

I don't  really  want  to  discuss  the  Correspondence  Theorem  or  its  applications. 
My  point  is  that  those  "memory  functions"  mentioned  in  it,  are  precisely  the 
equivalent  or  the  analogue  of  those  "complex-elastic  constants"  that  some  authors 
try  to  introduce  fictitiously  in  real-space  rather  than  in  the  Laplace-space , be- 
fore inversion  to  the  real  time-domain.  Hence,  the  way  to  use  those  "complex 
elastic  constants"  in  a way  that  H/orks,"  is  really  in  the  light  of  the  correspond- 
ence theorem.  Unfortunately,  that  is  quite  a difficult  task. 

We  have  already  pointed  out  that  another  method  which  "works"  is  to  measure 
the  (real)  elastic  and  viscous  constants  of  the  material,  each  set  independently 
of  the  other,  and  then  solve  the  field-equations  with  the  damping  terms  included, 
with  those  numerical  values  measured  for  the  constants.  This  method  is  also 
general  and  works  for  any  viscoelastic  model  and  it  is,  incidentally,  the  way 
viscous  flow  problems  are  attacked  in  Fluid  Mechanics. 


7 See  A.  Ceaal  Eringen,  "Continuum  Mechanics,"  John  Wiley  and  Sons,  Inc.,  1967, 
Article  9.12,  p.  368. 
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Observation  of  Table  1 and  2 shows  why  the  telegraph-type  equations  we  have 
discussed  here  are  so  important  in  viscoelasto-dynamieal  problems.  Note  that  the 
telegraph-type  equations  for  the  potentials  ip,  x are  of  the  type  given  here  in 
eq.  (l),  for  the  Kelvin-Voigt  model  only.  Note,  for  exmaple,  that  for  the  Maxwell 
model  they  are  substantially  different.  (Table  1,  eqs.  C-10).  The  equation  for 
<}>  is,  ^ 


26  + 
» 


hcr-§- 


4-1  v2$  = ^-  ft 

3tJ  ==  it 


^4  + 26  (3a  + 26) 
3t 


3*1 

■Stl 


(11) 


which  is  of  the  same  "type"  but  not  quite  of  the  form  given  here  in  eq.  (l). 


We  can  conceptually  always  introduce  complex  moduli  (as  in  eqs.  (7)  or  (10) 
or  their  analogues  for  the  "standard"  model)  provided  we  work  in  Laplace  space  as 
prescribed  by  the  Correspondence  Theorem.  If  we  introduce  them  in  the  real  time- 
domain  and  then  try  to  physically  identify  the  real  (or  imaginary)  parts  of  these 
complex  constants  with  the  real  elastic  constants  (or  the  viscous  coefficients), 
then  the  process  only  "works"  in  the  Kelvin-Voigt  model.  In  this  sense,  "works" 
means  that  the  equivalence  of  both  approaches  can  be  established. 


One  way  to  experimentally  measure  the  absorption  losses  of  material  samples 
is  the  impedance  tube  or  pulse  tube.  It  is  customary  for  the  literature  on  this 
technique  to  report  measurements  of  quantities  such  as  E'  and  E"  or  y'  and  y"  etc. 
This  technique  is  a good  source  for  the  popularity  of  "complex  moduli".  We  should 
note  that  the  lossy  samples  tested  in  this  fashion  are  always  characterized  by  the 
Kelvin-Voigt  viscoelastic  model  with  harmonic  time-dependence,  an  assumption  that 
may  not  always  be  Justified. 

2.  Conclusions 


We  summarise  our  points  as  follows, 

(a)  The  continuum  (infinite  number  of  degrees-of-freedom)  approach  to  viseo 
elasticity  as  a field-theory  is  the  only  way  to  go  today.  Discrete  approaches 
leave  much  to  be  desired  and  are  simplistically  unrealistic. 

(b)  It  looks  like  for  some  time  to  come,  we  will  be  dealing  with  the  Kelvin 
Voigt  and  Maxwell  models,  since  any  other  model  presents  too  many  analytical 
difficulties. 

(c)  A systematic  way  to  set  up  viscoelastic  boundary-value-problems  and 
solve  the  field  equations  in  these  two  models,  is  presented  here  for  cylinder  and 
sphere  problems.  This  approach  Is  shovn  in  Tables  1 and  2 which  are  self- 
explanatory. 

(d)  The  most  general  method  available  to  soi/e  static  or  dynamic  problems 
of  viscoelasticity  in  any  model  is  by  means  of  the  Correspondence  Theorem.  Here 
ve  must  work  in  Laplace-space  and  the  inversions  are  hard.  This  approach  is 
totally  equivalent  (at  least  for  the  Kelvin  and  Maxwell  models)  to  our  approach 
described  above  in  (3)  and  in  Tables  1 and  2. 
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(e)  Another  technique,  which  holds  for  any  model,  is  also  given  here.  It 
consists  of  solving  those  resulting  "telegraph-type"  equations  for  the  scalar 
potentials  keeping  in  them  the  damping  terms,  hut  using  real  values  for  the 
elastic  and  viscous  constants  which  are  to  be  experimentally  measured  as  it  is 
done  in  Fluid  Mechanics. 

(f)  We  claim  this  technique  is  more  systematic  and  less  confusing  than  to 
ignore  the  damping  terms  introducing  instead  "complex  elastic  moduli".  The 
complex  moduli  "trick"  only  works  for  the  Kelvin  model  with  harmonic  time- 
dependence  anyway. 

(g)  It  seems  clearer  to  us  when  talking  about  elastic  constants,  say, 

Young's  modulus  E,  to  think  of  the  slope  of  the  stress-strain  curve  as  found  in  a 
tension  test,  than  of  that  "complex  Young's  modulus"  which  contains  the  "loss"  in 
its  imaginary  part,  all  because  some  want  to  solve  wave  equations  rather  than 
telegraph  equations,  particularly  in  the  only  situation  when  one  is  no  harder  to 
solve  than  the  other. 

(h)  The  determination  of  the  solid's  elastic  constants  should  be  kept  sep- 
arate and  independent  from  the  determination  of  the  viscosity  coefficients.  I 
know  this  can  be  done  in  some  instances,  but  I am  not  aware  of  how  plausible  this 
recommendation  can  be  in  all  instances. 

(i)  Finally  it  is  worth  stating  that  the  pulse  tube  measurements  implicitly 
describe  the  viscoelastic  "losses"  in  the  sample  by  the  Kelvin-Voigt  model,  an 
observation  that  escaped  me  (and  perhaps  others)  until  recently. 

(J)  Equation  (10)  can  be  used  to  express  the  complex  shear  and  dilational 
moduli  for  the  Maxwell  model  in  terms  of  the  viscosity  coefficients  and  the  elastic 
constants  of  the  model  as  follows. 


These  relations  may  be  useful  when  trying  to  interpret  "losses"  in  the  light 
of  the  Maxwell  model  if  one  day  one  wishes,  or  the  need  arises  to  do  so.  It  Is 
obvious  that  these  relations  are  considerably  sere  complicated  than  the  analogous 
ones  for  the  Kelvin  model,  which  are, 

M*  a U £l  * IujH]  , .\*  ♦ Z\i*  =(A  +2u  ) [i  +i(iiM]  where, 

6 6 C 
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